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Abstract 

We present a complete calculation of the photon + 1 jet rate in e + e~ annihila- 
tion up to 0{aa s ). Although formally of next-to- leading order in perturbation theory, 
this calculation contains several ingredients appropriate to a next-to-next-to-leading 
order calculation of jet observables. In particular, we describe a generalization of the 
commonly used phase space slicing method to isolate the singularities present when 
more than one particle is unresolved. Within this approach, we analytically evalu- 
ate the singularities associated with the following double unresolved regions; triple 
collinear, soft/collinear and double single collinear configurations as well as those from 
the collinear limit of virtual graphs. By comparing the results of our calculation 
with the existing data on the photon + 1 jet rate from the ALEPH Collaboration 
at CERN, we make a next-to-leading order determination of the process-independent 
non-perturbative quark-to-photon fragmentation function D q ^(z, up) at 0(aa s ). As 
a first application of this measurement allied with our improved perturbative calcula- 
tion, we determine the dependence of the isolated photon + 1 jet cross section in a 
democratic clustering approach on the jet resolution parameter y cut at next-to- leading 
order. The next-to-leading order corrections to this observable are moderate but im- 
prove the agreement between theoretical prediction and experimental data. 



1 Introduction 



Although the production dynamics of photons and of jets are both individually well under- 
stood theoretically, their interplay is not. Hadronic jets are initiated by the production of 
quarks and gluons which subsequently fragment into clusters of hadrons. In events where a 
photon is produced in addition to the jets, this photon can have two possible origins: the 
direct radiation of a photon off a primary quark and the fragmentation of a hadronic jet 
into a photon carrying a sizeable fraction of the jet energy. While the former direct process 
can be calculated within the framework of the Standard Model, the latter is described by 
the process-independent quark-to-photon and gluon-to-photon fragmentation functions [3J, 
which cannot be calculated using perturbative methods but have to be determined from 
experimental data. 

Directly emitted photons are usually well separated from all hadron jets produced in a 
particular event, while photons originating from fragmentation processes are primarily to be 
found inside hadronic jets. Consequently, by imposing some isolation criterion on the photon, 
one is in principle able to suppress (but not to eliminate) the fragmentation contribution to 
final state photon cross sections, and thus to define isolated photons. However, recent analy- 
ses of the production of isolated photons in electron-positron and proton-antiproton collisions 
have shown that the application of a geometrical isolation cone surrounding the photon does 
not lead to a reasonable agreement between theoretical prediction and experimental data. 
This discrepancy has important consequences. Data on high photon production are com- 
monly used to extract informations on the gluon distribution in the proton 0. Furthermore, 
the understanding of final state photon radiation in proton-proton collisions is even more 
crucial for searches for the two photon decay mode of the Higgs boson at the LHC ||. 

An alternative approach to study final state photons produced in a hadronic environment 
is obtained by applying the so-called democratic clustering procedure ||. In this approach, 
the photon is treated like any other hadron and is clustered simultaneously with the other 
hadrons into jets. Consequently, one of the jets in the final state contains a photon and is 
labelled 'photon jet' if the fraction of electromagnetic energy within the jet is sufficiently 
large, 

Eem ^ ,.. 1 x 

Z = p — p > Z cut, (1-1) 

&EM + &HAD 

with z cnt determined by the experimental conditions. This photon is called isolated if it car- 
ries more than a certain fraction, typically 95%, of the jet energy and said to be non-isolated 
otherwise. Note that this separation is made by studying the experimental z distribution 
and is usually such that hadronisation effects, which tend to reduce z, are minimized. The 
cross section for the production of isolated photons defined in this approach thus receives 
sizeable contributions from both direct photon and fragmentation processes. 
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Up to now, this democratic procedure has only been applied by the ALEPH collaboration 
at CERN in an analysis of two jet events in electron-positron annihilation in which one of the 
jets contains a highly energetic photon ||. In this analysis, ALEPH made a leading order 
determination of the quark-to-photon fragmentation function by comparing the photon + 1 
jet rate calculated up to 0(a) @] with the data. Then, by inserting their measurement 
into an 0(aot s ) calculation of the 'isolated' photon + 2 and 3 jet rates, they found good 
agreement for a wide range of values of the jet resolution parameter y cut . 

In addition, the 0(a) 'isolated' photon + 1 jet rate agreed well with the ALEPH data. 
This is remarkable, since with the use of an isolation criterion based on the application of 
a geometrical cone previously considered ||, the agreement between data and theory for 
this rate was poor and only improved at the expense of large negative radiative corrections 
[0, ||, [| |10| . The main reason for the size and sign of these corrections was due to the 
exclusion of soft gluons from the cone containing the photon |]TT . 



Recently, we have presented a next-to-leading order determination of the quark to pho- 
ton fragmentation function ]12]| based on an 0(aa s ) calculation of the photon + 1 jet rate 



and the ALEPH data. Both theoretical and experimental analyses were performed using 
a democratic clustering approach; in particular both analyses used the Durham jet recom- 
bination algorithm [13[]. Not only is the agreement with the 'isolated' photon + 1 jet rate 



improved, but the higher order corrections obtained in this approach were found to be mod- 
erate, demonstrating the perturbative stability of this particular photon definition. In this 
paper, we describe the calculation in more detail. In particular, we extend the commonly 
used phase space slicing technique [14], [L^] developed to isolate the divergences when one of 



the final state particles is unresolved to situations where two particles are unresolved. 

The paper is organized as follows. In section [| we discuss the general structure of photon 
+ n jet events in electron-positron collisions, with particular emphasis on where the singular- 
ities arise and how they are absorbed into the photon fragmentation functions. In this fixed 
order approach, it is critical that the power counting of the various contributions is correct. 
One consequence is that the commonly held view |16| that the quark-to-photon fragmenta- 
tion function is 0(a/a s ) is untenable. The particular features of the photon + 1 jet rate 
are discussed in section |3|. We adopt the resolved parton philosophy of [ll| and introduce 
a small parameter y min to separate the final state phase space into resolved and unresolved 
regions and to isolate the singularities lying in the unresolved regions. Within this approach, 
it is vital that the various singular regions precisely match onto each other. The evaluation 
of the associated different divergent contributions to the photon + 1 jet rate is presented in 
sections Fully resolved and single unresolved contributions from the four parton process 
7* — > qq^g are described in section [| while the double unresolved contributions that arise 
when, for example, three particles are simultaneously collinear are analytically evaluated in 
section Contributions from the one-loop process with collinear final state partons are pre- 
sented in section |] and the fragmentation contribution is presented in section [7]. All divergent 
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contributions are combined in section || where we explicitly show that after factorization 
of collinear singularities, no divergences remain. Although the individual contributions are 
finite, they do depend on the parameter y m i n . Clearly, the physical cross section cannot 
depend on y m i n , and in section [| we numerically show that, provided y m i n is taken small 
enough, this is indeed the case. We combine our numerical results with the experimental 
data in section |10] and determine the next-to-leading order quark-to-photon fragmentation 
function for two values of a s [Mz)- Finally, section 11 contains our conclusions. 



2 The n jet + photon cross section 

Let us first consider the general structure of the e + e~ — > n jet + photon cross section, fully 
differential in all quantities, 

da(n jets + "7") = da y + £da a <g> D^ r (2.1) 

a 

There are two contributions, first the 'prompt' photon production where the photon is pro- 
duced directly in the hard interaction and second, the longer distance fragmentation process 
where one of the partons fragments into a photon and transfers a fraction of the parent 
momentum to the photon. Each type of parton, a, contributes according to the process 
independent parton-to-photon fragmentation functions and the sum runs over all par- 

tons. Note that although the fragmentation functions are non-perturbative, we can nominally 
assign a power of coupling constants, based on counting the couplings necessary to radiate 
a photon. Since the photon couples directly to the quark, D q ^ is naively of O(o) while 
the gluon can only couple to the photon via a quark and we might expect that D g ^ is of 
0(aa s ). This simplistic argument is supported by models of the fragmentation function [|17|] 
which suggest that gluon fragmentation is a much smaller effect than quark fragmentation. 

The individual terms in eq. (|2.1| ) may be divergent and are denoted by hatted quantities. 
However, we can reorganise them as follows. In the most general case we have, 

d(T 7 = dcx 7 + ^2 & a q ® ^q^i + & a g ® r g ^ 7 , 
q 

da q = do Y ® ^q'->q + dog <g> T g ^ q , 
q' 

dcr g = d °g ® T q->g + dcr s ® r 9^<?> ( 2 - 2 ) 
q 

where all of the divergences are concentrated in the factorization scale dependent transition 
functions IV^-. The sum covers all active quark and antiquark flavours. The process specific 
cross sections da 7 , da q , da g are now finite. Inserting these definitions back into eq. ( ]2.1p , 
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yields a physical cross section in terms of finite quantities, 

da(n jets + "7" ) = da 7 + da q ® D q-*i + dcr 9 ® A?-7 ( 2 - 3 ) 

<? 

where the physical (and factorization scale dependent) fragmentation functions are given by, 

= r„ + Er„^ 9 l 7 + r„®D s l r (2.4) 



While the singular parts of the transition functions are process independent and well 
known, it is sometimes convenient to carry some finite perturbative pieces as well so 
that, 

da R = d<j 7 - ^2 da q <g> a q ^ - da g ® a g _> 7 , 

r R = r + o 

^9-^7 = ^9^7 + a 9-^7' (2-5) 

and to define the 'effective' parton-to-photon fragmentation functions, 

— >7 — > 7 ~t~ (2, 
^9^7 = "^9^7 "I" °9— 



£> 9 ^ 7 = -D 5 ^ 7 + a 5 _ 7 . (2.6) 



With these modifications that are particularly suited to the resolved parton philosophy [[L5[] 
we shall employ later on to isolate the singularities in d<r 7 and da a , the cross section is given 

by ' _ 

da(n jets + "7" ) = da* + ® ^9-7 + d(T f ® ^9-7 > ( 2 -7) 



where the superscript R indicates that the partons are 'resolved' |15| . Within the fixed order 
approach, only the 'prompt' photon process contributes at lowest order, and, 

da^ = daf LO) =edao(np + 7), 
d<xf°) = 0, 

daf^ = 0. (2.8) 

The n parton + photon cross section d<5"o(n p + 7) is evaluated in the tree approximation and 
G represents the experimental jet and photon definition cuts. In this way the theoretical 
cross section can be matched onto the specific experimental details. At this order each parton 
is identified jet and the photon cLS cL photon. 
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At next-to- leading order, both 'prompt' photon production and the quark-to-photon 
fragmentation process contribute, 

da (NLO) = @ P + l) + J da ((n + 1) p + 7)) , 

da (NLO) = Q d ( 7o((n + l)p), 

d& {NLO) = q (29) 

'Prompt' photon production may occur via the one loop virtual corrections to the n parton 
+ photon process d<3"i(n p + 7) or by the tree level emission of an additional parton d&o((n + 
1) p+7). The integral sign represents the integration over the additional phase space variables 
due to the presence of an extra parton in the final state. The fragmentation contribution 
arises from the lowest order n + 1 parton process with one quark fragmenting into a photon. 
The remaining unfragmented partons are directly identified as jets. Although the physical 
cross section is finite, individual contributions are divergent. In da^ NLO \ there are infrared 
singularities arising from configurations where one of the partons is theoretically 'unresolved'. 
For example, the virtual graphs contain singularities due to soft gluons or collinear partons 
while similar divergences occur in the bremsstrahlung process. The correct treatment of 
infrared divergences is well known |Tj| [T!| and has been discussed widely in the literature. 
Many general approaches have been developed to isolate the infrared divergences at next-to- 



leading order [15], g(| £y, |22fl. In J|, the approach of |15j was used to isolate the divergences 
present in the bremsstrahlung contributions and to combine them with those arising from 
the virtual graphs. While the leading infrared singularities due to soft gluon emission should 
cancel within da^ NLO \ some collinear mass singularities remain. These singularities, which 
occur as the quark and photon become collinear, are proportional to da^ NLO \ are factorizable 
and can be absorbed by a redefinition of the fragmentation function as described above. 
Working within the MS scheme ||23|| , the transition functions for a quark of charge e q are 
up to 0(a) given by, 



r„ 




'47T/Z 2 



e 3^7 



2tt ; r(i-e) v /4 J L e' q ^\ 1 v2 ' 10) 

where -P„^L is the e — > part of the the lowest order splitting function in (4 — 2e)-dimensions 



H, 

p„w = 1 + (1 -;' 2 - 2 . (2.1D 

Here, the variable z represents the fraction of the 'photon' energy that is genuinely electro- 
magnetic in origin in the quark-photon cluster. 

Each term in the 'prompt' photon contribution is finite and may be evaluated numerically, 
da R(NLO) = q p + 1 ) + J da R((n + i) p + 1 )\ (2.12) 
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where the superscript R indicates that the partons are 'resolved' ||15|| . For the numerical 
evaluation of the cross section within the resolved parton approach it is moreover useful 
to identify a universal constant piece and to absorb it into the 'effective' fragmentation 
function Q, 

2W.) = IW*, A*) + (f?) (log ( S ^*w%~ Z) ) P ^ (Z) + Z ) ' (2 ' 13) 

where M is the mass of the final state and s qq - is the invariant mass of the radiating quark- 
antiquark system. We see that the effective quark-to-photon fragmentation function, T> q ^ 7 , 
depends on the unphysical parton resolution parameter s min introduced to isolate the sin- 
gularities. Physical cross sections cannot depend on s min , and there must be a cancellation 
between the resolved prompt photon and fragmentation contributions. In ref. [|J this was 
shown to be the case. The genuine non-perturbative fragmentation function D q ^ 7 (z, [ip) 
cannot be calculated and needs to be extracted from data. 

At next-to-next-to- leading order, the 'prompt' photon contribution now contains contri- 
butions from the two-loop n parton + photon process, as well as the one-loop (n+ l)-parton 
+ photon and tree level (n + 2)-parton + photon processes. Similarly, the quark-to-photon 
fragmentation term receives contributions from the one-loop (n + l)-parton process and the 
tree level (n + 2)-parton process, while the gluon-to-photon fragmentation term, appearing 
for the first time at this order, only receives a contribution from the tree level (n+ l)-parton 
process, 

da (NNLO) = e^ 2 (7ip + 7)+y'da 1 ((7i + l)p + 7) 

+ J J d«7o((n + 2) V + l) 
d& (NNLO) = Q ^((n + i) p ) + y da ((n + 2) p)) , 
da (NNLO) = edoo((n+l) p). (2.14) 

As at next-to- leading order, the singularities present in the 'prompt' photon bremsstrahlung 
and in the fragmentation processes can be isolated. However, unlike at next-to-leading 
order, there are contributions to the 'prompt' photon bremsstrahlung process where two of 
the final state partons can be simultaneously theoretically unresolved. These contributions 
are expected to appear in the divergent part of any jet cross section evaluated at next- 
to-next-to-leading order. In the most general case, the following four double unresolved 
configurations can occur: 

• The triple collinear configuration, where three partons become simultaneously collinear 
to each other, but none are soft. 
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• The soft/collinear configuration, where two partons become collinear while a third 
parton becomes soft. 

• The double single collinear configuration, where two distinct, independent pairs of 
partons become simultaneously collinear. 

• The double soft configuration, where two partons become simultaneously soft. 



The photon + 1 jet cross section at 0(aa s ), calculated in the remainder of this paper, 
receives only contributions from three of these double unresolved configurations (triple 
collinear quark-photon-gluon, soft gluon/collinear photon and double single collinear quark- 
photon/antiquark-gluon) and the corresponding double unresolved factors will be calculated 
analytically in section [5[ There are also singular contributions arising when a gluon is vir- 
tual while another parton is collinear or soft. In our calculation, we have only a virtual 
collinear contribution, studied in section ||], as the photon cannot be soft. Finally, another 
new class of singularities occurs when together with the remnants of the quark/photon frag- 
mentation cluster an unresolved or virtual gluon is emitted. We will explicitly evaluate this 
fragmentation contribution in section [7]. Having isolated the singular contributions in this 
way allows us to define finite 'resolved' parton cross sections for both 'prompt' photon and 
fragmentation processes which in the most general case reads, 

d(J R(NNLO) = q ^ d(J R( n p + J da R ((„ + 1) p + 7 ) 

+ JI da*((n + 2)p + 1 ) 

da R(NNLO) = Q Ua?{{n + 1) p) + J da*((n + 2) p)) , 

da R(NNLO) = edaf((n + l)p). (2.15) 

Each resolved cross section is finite and can be numerically evaluated for arbitrary experi- 
mental constraints. 



The most singular divergences are due to soft gluon radiation and cancel between the 
virtual and real processes. However, double and single poles in e due to collinear singularities 
remain, whose form is determined by the known next-to-leading order transition functions. 
Keeping only terms that are necessary to expand eq. ( |2.2|) up to 0(aa s ) and bearing in mind 
that D g ^ is already 0(aa s ), then in the MS'-scheme, 



<?^7 



2tt ) r i 




T 2 (l 



.Ip(o) 
{ /4 , 



2e r 



J_P(0) 
2e 2 1^1 



P(0) _ 1_ p(l) ' 

2e q ~* 1 



7 



ae l q \ fa s T R 



g^i 



g^g 



2tt J V 2tt J T 2 (l 



a s C F 



. - e) V fi 2 F J W 2Z t 2e 7 . 



1 + 



1 /47T/1 2 

2rr 7 r(l-e) ^^4" 

1 /47T/i 2 V 



.Ip(o) 



2tt y r(l-e) V /4 



- p(O) 



0. 
1. 



(2.16) 



where the QCD casimirs C F and Tr are, 

iV 2 - 1 

C F -- 



2N 



T, 



R 



The lowest order quark-to-quark and gluon-to-quark splitting functions are [^4 

1 + z 2 3 



(1-*)+ 2 
, 2 + (l-,) 2 . 



5(1 - z) 



(2.17) 
(2.18) 



while the next-to-leading order quark-to-photon and gluon-to-photon splitting functions are 
given by ^ f2G|, 

1 9 
2 



+ 



ln 2 (l - z) + 41nzln(l - z) +8Li 2 (l - z) - -n 2 



P^%(z), (2.19) 



82 46 / 8 8 \ 

P$L{z) = -2 + Qz-—z 2 + — +[5 + 7z + -z 2 + —)\nz+(l + z)\n 2 z. (2.20) 

9 9z \ 6 ozj 

The transition functions determine the singularity structure of the bare fragmentation 
functions up to 0(aa s ), 



D B (z) 



<?^7 



Dg^ y (Z,jJ F ) + 



ae: 



1 /47T/i 2 \ 1 



_p(0) 



+ 



+ 



a s C F 



2tt / T(l - 
'ae 2 \ fa s C F 



2tt J r(l - e) V /4 / e 



1 / 47T/i 

2tt y v 2tt y r 2 (i - e ) \ i4 



2\ 2e r 



1 



1 



_p(0) (o, p(0) j p(l) 

2 e 2 ^ 2e 
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(2.21) 



Dg^y{Z,(J, F ) 



E 



2vr 



2tt / r 2 fl - el 



'4:71/1 



2\ 2e r 



— P (0) (5?) P (0) 4- — P (1) 
2 e 2 s-* 1 ? ?~*7 2e 9_>7 



'2.22) 



which must absorb the remaining collinear singularities present in the 7 + n jet cross section 
at next-to-next-to-leading order. Note that in expanding the perturbative counterterms, we 
have systematically assumed that the gluon-to-photon fragmentation function is 0(aa s ), 
and therefore does not contribute to Dg at this order. 

Concentrating on the evaluation of the photon + 1 jet rate at fixed order, 0(aa s ), the 
production of gluons is suppressed by a power of a s compared to the quark production. 
Therefore, the contribution from the gluon fragmentation to the photon + 1 jet rate is of 
0(ao> 2 s ) and must be neglected in a consistent fixed order framework. Consequently, the 
quark-to-photon fragmentation function alone must absorb the remaining singularities and 
demanding that this is so provides a powerful check that the singularities have been correctly 
isolated. This factorization procedure will be explicitly described in section [8]. Moreover, 
we note that evaluating the real double unresolved contributions using a strong ordered 
approach, i.e. requiring that one unresolved parton is much less unresolved than the other, 



fails to produce the necessary singularity structure P? . 



3 The 'photon' + 1 jet cross section 

We now concentrate on the case where only a single jet is produced in addition to the 
'photon'. The photon + 1 jet rate in electron-positron annihilation is rather particular, 
since the leading order process, 1 parton + 'direct' photon production, vanishes. The lowest 
non- vanishing order ||, 0(a), is therefore equivalent to the next-to-leading order outlined 
above, where 'direct' and fragmentation processes contribute at equal level. This observable 
is therefore particularly suitable for determining the non-perturbative component of the 
quark-to-photon fragmentation function. 

The parton level subprocesses contributing to the photon + 1 jet rate at 0(aa s ) are 
shown in Fig. |l[ 

(a) The tree level process 7* — > qqg'j, where the final state particles are clustered together 
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Figure 1: Parton level subprocesses contributing to the photon + 1 jet rate at 0(aa s ). 



such that a "photon jet" and one additional jet are observed in the final state. As the 
photon must be identified in the final state it cannot be soft. 

(b) The one loop gluon correction to the 7* — > qqj process, where the photon and one of 
the quarks are clustered together or the photon is isolated while both quarks form a 
single jet. 

(c) The process 7* — > qqg, where one of the quarks fragments into a photon while the 
remaining partons form only a single jet. 

(d) The one loop gluon correction to the 7* — > qq process, where one of the quarks frag- 
ments into a photon. 

(e) The tree level process 7* — > qq with a generic 0(aa s ) counterterm present in the bare 
quark-to-photon fragmentation function. Inclusion of this contribution absorbs all left 
over singularities of the processes cited above. 



Each of these processes gives a singular contribution to the final photon + 1 jet rate. 
Therefore, to combine them together numerically in a way that can match onto the precise 
experimental cuts, we must first isolate the divergences analytically, so that the remaining 
finite contributions may be numerically evaluated. We choose to employ the hybrid sub- 
traction method [2S| which extends the phase space slicing method described in []r5| |. In 
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this approach, we introduce a parton resolution parameter y m [ n to decide whether or not two 
partons are theoretically resolved; if, for example, < y m i n , then partons i and j are not re- 
solved. In these unresolved regions, the matrix elements are singular, typically proportional 
to 1/yij, and we analytically integrate approximate forms for the exact matrix elements over 
the restricted phase space. The essence of this approach can best be illustrated using the 



simple one-dimensional integral as suggested by Kunszt and Soper in |29 





I = ]im{ —x e F(x)--F(0)\. (3.i; 



F(x) is a complicated function, which renders the evaluation of X analytically impossi- 
ble. X could represent a n-jet cross section while F(x) could stand for the (n + l)-parton 
bremsstrahlung matrix elements and x for a scaled invariant mass y^. As x — > 0, which 
corresponds to the case when one of the final state particles becomes soft or collinear, the 
integrand is regularized by the x e factor as in dimensional regularization. The first term is 
however still divergent as e — > 0. This divergence is cancelled by the second term - which is 
the equivalent of the n-parton one-loop contribution - so that the integral is finite. 



In the hybrid subtraction method |28| , we choose to evaluate the integral by adding and 
subtracting the approximate matrix elements denoted by F(0) in the unresolved regions of 
phase space, so that, 

{r 1 dx rVmin dx rVxmn dx I 1 

/ ™ x 'F(x) - F(0) / — x e + F(0) / — x e - -F(0) . 
Jo x Jo x Jo x e J 

The approximate matrix elements are both simpler and process independent and the inte- 
grations over the unresolved phase space can be carried out analytically, 

2/min dx , 1 



so that, 



/•3/min CLX 

F(0) / -x* = ~F <4 in , 
Jo x e 

l im / f 1 — X -F(x) + T min — x e \F(x) - F(0)1 
^0 {Jy min x Jo x 

+lF(0)^ in -iF(0)} 



f 1 dx , r Vmin dx 

'Umin X 



F( x ) + f Vmin ^ [F(x) - F(0)] + F(0) \n(y min ) } . (3.2) 
jo x 



All three terms are now finite so that the e — > limit may be safely taken and, furthermore, 
are in a form suitable for numerical evaluation. 

Of course, the parton resolution parameter y min is unphysical and the integral X or 
equivalently any physical process should not depend on y m i n . The y min dependence of the 
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three terms in eq. ( |3.2|) should therefore cancel. In the evaluation of a jet cross section, 
this cancellation is usually realised numerically by a Monte Carlo program. This is not a 
straightforward point for the following reasons. The matrix element approximations used 
in the analytic part of the calculation are reliable only when y m i n is small and are best 
when y m j n is the smallest possible. On the other hand, smaller values of y m i n generate larger 
cancellations amongst the terms, possibly giving rise to numerical instability problems. In 
practise y m i n is chosen in such a way that the approximations performed in the analytic 
calculation are valid and that the numerical errors are minimized. 

This approach, or the more basic phase space slicing method where the second term 
in eq. ( ^.2| ) is neglected, has been applied to a wide variety of processes; e + e~ — > 2 jets, 
e + e~ — > 3 jets |15| , pp — > W, Z + 1 jet, pp — > 2 jets |3(J and ep — > e + 2 jets f3~If . In all of 



these next-to-leading order calculations, there can be at most one unresolved particle; one 
soft gluon or two collinear partons. In our calculation of the photon + 1 jet rate at 0(aa s ), 
we will be concerned, for the first time, with situations where two partons are unresolved. 

So far we have described the various parton level processes contributing to the photon + 1 
jet rate at 0(aa s ) and merely outlined the general structure of the contributions associated 
with these processes. Each of these processes contains different singular contributions that 
arise when one or more particles are theoretically unresolved. We now consider each process 
in turn, and use the parton resolution parameter y min to define for each of them, the resolved 
and unresolved contributions from the resolved and unresolved regions of the phase space 
respectively. As we will see in sections all theoretically unresolved contributions have 
the common structure already encountered in fll5f . They can be written as the product of 



an unresolved factor containing all the singularities and a resolved cross section. 

Note that in isolating the singularities, we are not concerned with the precise experimental 
jet definition. Provided the theoretical resolution parameter y m i n is sufficiently small, any 
jet algorithm can be imposed numerically on the resolved parton cross sections, da R . 



3.1 7* — > qq-j with real gluon bremsstrahlung 

The various configurations where the tree level process 7* — > qqg^ contributes to the photon 
+ 1 jet rate are illustrated in Fig. [| Note that topologies where the role of quark and anti- 
quark are exchanged are also present, but are not shown. Jets formed when the experimental 
jet algorithm clusters particles together are denoted by (ij), while theoretically unresolved 
clusters are denoted by [ij]. The associated real contributions can be separated into three 
categories: either theoretically resolved, single unresolved or double unresolved. Care must 
be taken to divide the phase space so that the different regions match onto each other. We 
must ensure that the difference of exact and approximate matrix elements tends to zero as 
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(c) 
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([97]) (t9gl) 


(i) 
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Figure 2: Different final state 'photon' + 1 jet topologies arising from the tree level 7* — ■> gf^yg 
process. The 'photon' jet is moving to the left while the recoiling hadronic jet moves to the 
right. Square brackets denote theoretically unresolved particles, round brackets experimental 
clusters. 



the singularity is approached and equally, there is no overcounting of the singularities. This 
can be done by constraining the variables, 

2/975 2/975 2/995 2/995 
although in the double unresolved regions, we shall choose to constrain the combinations, 

2/9795 2/9795 

for certain configurations. The matrix elements contain poles in each of these six invariants, 
and are therefore singular as any of them tends to zero. There are two invariants, y gi and y q q, 
which are not associated with any singularities and are therefore completely unconstrained. 

The individual configurations can be structured as follows: 
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Theoretically resolved contributions 

If all particles are resolved, a 7 + 1 jet event can only be formed if some final state par- 
ticles are clustered together by the jet algorithm. The possible configurations yielding 
a photon + 1 jet event are displayed in Fig. |2|.a. In principle this region is defined by re- 
quiring that all invariants are greater than the theoretical parton resolution parameter, 
y min , i.e. by requiring that, 

y<Vi ^* Vm'mi V(T1 ^* Vmmi Vqg Vminj Vqg Z/miri' (^"^) 

However, it turns out that the boundaries of this region are more subtle than that and 
must be chosen to match onto the boundaries of the unresolved regions. Consequently, 
the fully resolved region is defined as being the remaining phase space region of the 
four parton phase space when all unresolved regions are excluded. 

Single unresolved contributions 

As shown in Fig. §.b, there are two classes of single (or one-particle) theoretically 
unresolved contributions, depending whether the gluon or the photon is unresolved. 
The single unresolved regions associated to these contributions are defined as follows: 

(i) The collinear quark-photon region 

If the photon is unresolved, it is collinear to the quark while the gluon is hard, 
i.e. the gluon is theoretically resolved but combined with the photon-quark cluster 
or with the antiquark by the experimental jet algorithm. Alternatively, the gluon 
forms a jet on its own while the antiquark is clustered into the photon jet. This 
can be defined by the constraints, 

Vq~1 ^ Vmm Vqq-yi Uq-y > VminUqq')-, Uq'yg ^* Urn'mi Vqg > 2/min; 

(3.4) 

where y qq -j is the scaled invariant mass of the radiating quark-antiquark-photon 
antenna, s q q^/M 2 . 

(ii) The unresolved gluon region 

If the gluon is theoretically unresolved, it can be soft or collinear to the quark or 
to the antiquark, while the photon is experimentally combined with the quark to 
form the photon jet or is isolated while all other partons form a single jet. The 
three possible regions are defined by: 

(1) The collinear quark-gluon region 

Vqg < Urn'm Vqqg i Vqg > UmmU qqg: Uq'yg ^* Vminj Vqj > 2/min; (3.5) 

(2) The collinear antiquark- gluon region 

Vqg < Urn'm Vqqg i Vqg > Urn'm Uqqg, Uq'yg ^* Ummj Vqj > 2/min; (3.6) 
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(3) The soft gluon region 



Uqg ^ yminUqqg-i Uqg ^ ?/min Dqqg > Vqf ^* 2/mim 2/<j>7 ^* Z/min- 

As in the unresolved photon case, the mass of the antenna radiating the gluon, 
y q qg plays a role in determining when the approximate matrix elements are used. 

(c) Double unresolved contributions 

These contributions arise when both the photon and the gluon are theoretically "un- 
seen" in the final state. As the photon has to be seen in the final state, it can only 
be collinear with the quark and cannot be soft. The gluon on the other hand can be 
collinear with the quark or with the antiquark or it can be soft. Corresponding to 
these different final state configurations we define three double unresolved phase space 
regions: 

(i) The triple collinear region 

The photon and the gluon are simultaneously collinear with the quark. Here we 
need to constrain the "triple" invariant y qjg = y qi + y qg + y ig since it appears in 
the denominator of the four-particle matrix elements, 

and yq g > y min . (3.8) 

Note that this implies y qi < y min and y qg < y m i n . We require y qg > y min since in 
this region the gluon is collinear but not soft. 

(ii) The soft/ collinear region 

The photon is collinear with the quark while the gluon is soft, 

and y qg < y min and y qg < y min . (3.9) 

(iii) The double single collinear region 

The photon is collinear with the quark while the gluon is collinear with the anti- 
quark, 

and y qg > y min and y qg < y min . (3.10) 

For these three contributions, the final state configuration corresponds already to a 
photon +1 jet event. Hence, the final state particles will not be clustered further by 
the jet algorithm. These contributions are schematically displayed in Fig. |2].c. As 
before, these regions are matched by three analogous double unresolved regions where 
the photon clusters with the antiquark. 

The decomposition of the four-particle phase space is summarized in Fig. |3|. In this table, 
we have specified which invariants are less than y min (or a cut proportional to y m m) for each 
singular region of phase space. We have also noted which invariants are greater than y m[n to 
eliminate overlaps between regions determined by the same combinations of invariants less 
than y m { n . Invariants that are not specified are completely unconstrained. 
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Real contributions 



Resolved 



Single unresolved 



q'-f collinear 



Double unresolved 



qg collinear 



qg collinear 



g soft 



Vq-f ^ 2/min2/997 
Vq-fg ^* 2/min 
Uqg ^ 2/min 

2/g~7 ^* 2/min2/997 



Vqjg ^* 2/min 
2/95 ^ 2/min2/ggg 
Uqg ^ 2/min2/ggg 
Uq-/ ^* 2/min 



2/<JT ^ 2/min 
2/g<? ^* 2/min2/ggcj 
2/gg 2/min2/ggg 
Vq.19 Drain 



Vq-f ^ 2/min 
2/gg ^ 2/min2/ggcj 
Uqg < ymmUqqg 

Vti - > 2/min 



Triple collinear 



Soft / collinear 



Double single collinear 



2/g79 2/min 
2/g~9 ^* 2/min 
2/g~7 ^ 2/min 



2/?7 ^ 2/min 
2/(79 ^ 2/min 
2/<79 2/min 
2/<?7 ^ 2/min 



2/?7 ^ 2/min 
2/(79 ' > 2/min 
2/99 ^ 2/min 
2/<?7 ' > 2/min 



Figure 3: Phase space decomposition of the real 7* — > qqjg contributions. Note that the 
single and double unresolved regions where the photon clusters with the antiquark are not 
shown. For these regions, the necessary cuts are obtained by exchanging q and q. Altogether, 
there are five single unresolved and six double unresolved regions. 



3.2 7* — > qqj with a virtual gluon 

The final state topology is similar to that for the lowest order 7* — > qq^ contribution dis- 
cussed in (U and contributes to the photon + 1 jet rate, if two of the final state partons 
are clustered together. We distinguish the theoretically unresolved collinear quark-photon 
contribution from the contributions where the theoretically resolved photon is clustered with 
the quark by the experimental jet algorithm to form the photon jet or isolated while quark 
and antiquark combine to form the jet. The three particle phase space therefore divides as 
follows: 



(a) The resolved photon region 
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(b) The quark-photon collinear region 

2/<J7 *^ 2/min; Vq-y ^* Z/mim (^-12) 

plus a similar region for the collinear antiquark-photon configuration. 



3.3 7* — > qq(g) with fragmentation 

The 0(a s ) processes with associated fragmentation shown in Fig. [Tj.c and |l|.d contribute to 
the photon + 1 jet cross section if, in addition to the photon-jet, only a single jet is formed. 
As illustrated in Fig. |], the gluon may be resolved, unresolved or virtual, while the photon 
is produced via the fragmentation process. There are five distinct contributions: 

(i) The resolved gluon region 

The real gluon is theoretically resolved, but may be clustered by the experimental jet 
algorithm with either the photon fragmentation cluster or with the antiquark or it may 
form a jet on its own, while the antiquark is combined with the photon/quark cluster. 

Uqg ^* 2/miru Uqg ^* 2/min- 

(3.13) 

(ii) The unresolved gluon region 

If the gluon is theoretically unresolved, it can be soft or collinear to the fragmenting 
quark or to the antiquark. The three possible regions are defined by: 

(1) The collinear quark-gluon region 

Uqg *^ 2/min; Uqg > Vmin- 

(2) The collinear antiquark- gluon region 

Uqg V 'min; Vqg ^ Vmin- 

(3) The soft gluon region 

Uqg *^ 2/min; Vqg ^ 2/min- 

(iii) The gluon is virtual. 

Since this process is already of 0(a s ), only the 0(a) counterterm in the bare fragmentation 
function contributes. 
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(3.14) 
(3.15) 
(3.16) 



(i) 


(y) (qg) (yg) (q) (yq) (g) 


(ii) 


X' ^ — — 

(y) dqg]) dig]) ® (y) ® 


(iii) 


(y) ® 



Figure 4: Different finaf state 'photon' + 1 jet topologies arising from the 7* — > qq(g) process 
with subsequent fragmentation of the quark into a photon. The 'photon' jet is moving to the 
left while the recoiling hadronic jet moves to the right. Square brackets denote theoretically 
unresolved particles, round brackets represent experimental clusters. 



3.4 7* — > qq with fragmentation 

In addition to the processes described above involving a real or virtual gluon, one has to 
consider a contribution to the photon + 1 jet rate from the generic 0(aa s ) counterterm 
present in the bare fragmentation function of eq. ( |2.21| ). 

4 Resolved and Single unresolved contributions 

In the next two sections, we will discuss the contributions to the 7 + 1 jet rate at 0(aa s ) 
relevant to the process 7* — > qq^y with real gluon bremsstrahlung. As shown in Fig. |], there 
are many different topologies contributing which can be divided according to the number of 
theoretically unresolved particles. In this section, we consider the cases where at most one 
particle is unresolved, while the double unresolved contributions are studied in section || 

The 7* — > qqjg process contributes to the 7 + 1 jet differential cross section if the final 
state configuration is such that only the photon jet and a single associated jet are observed. 
The generic contribution of this four particle final state process to the rate can be written 
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as, 



Here |A4 g g 7S | represents the squared matrix elements in ^-dimensions, while the Lorentz 
invariant phase space for the decay of an off-shell photon with p 2 * = M 2 can be written, 

J dPi d) (M;p q ,p^,p g ,pg) = 2 7-4e(27r)8-6e J dfid ~ 1 drid - 2 dQ d-sdy q q dy g7 dyqgdyq 7 dyq g dy jg 

x (-A A y 1/2 ~ € 5(y qq + y qj + y qg + y^ +y qg +y 19 - 1), (4.2) 

with, 

^4 = yqq y~/g ~\~ 2/cJ7 Z/qg + Z/gg Vq'y 2 ( T/ggT/g-y?/^?/^ ~\~ V 'q"/V 'q"/V 'qgV 'qg + 2/gq2/(jg2/7g2/<J7 ) • (4-3) 



The matrix elements are well known and contain terms up to 0(e 3 ). These additional 
contributions, which vanish in 4-dimensions, are important in determining the contributions 
from the unresolved regions which may contain singularities up to (9(l/e 3 ). 



4.1 The resolved contributions 

In the resolved phase space region which is the region of the four-particle phase space left 
over when all unresolved regions specified in Fig. |3| are excluded, the matrix element squared 
is finite. Thus we may take the e — > limit for both matrix elements and phase space, so 
that 

d<4g g = (^) (^pj 4(2vr) 4 ^ / \Mq^\ 2 dPi R \M;p q ,p q , Mg ), (4.4) 

( 

where the phase space dP 4 is restricted to this resolved region. For these resolved contri- 
bution, the integration is carried out numerically and the photon and jet definitions applied 
directly to the final state particles. As we are employing the hybrid subtraction method to 
calculate the photon +1 jet cross section, there are additional contributions which are to be 
numerically calculated. Those come from evaluating the difference between the exact and 
approximate squared matrix elements in the various unresolved regions of phase space. 



4.2 The single unresolved contributions 

There are two classes of single unresolved real contributions depending on whether the photon 
or the gluon is unresolved. If the photon is unresolved, it is collinear to the quark while if 
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the gluon is unresolved it can be collinear to the quark, collinear to the antiquark or it can 
be soft. The possible final state configurations of single unresolved contributions yielding 
a 7 +1 jet event are displayed in Fig. @.(b). In each unresolved region of the four-particle 
phase space we expect to be able to write the differential cross section as the product of a 



one particle unresolved factor and a three-particle cross section |T5|| . 



4.2.1 The unresolved gluon contribution 

The phase space region where the quark and the gluon are collinear is defined by eq. Q3.5p, 



Uqg ^ Z/min Dqqg i Dqg ~> Drain Dqqg i Vqig V 'min; Vq-y ^* 2/min- 

Here, the quark and the gluon cluster to form a new parton Q such that, 

Pq+Pg =PQ, 

and carry respectively a fraction y and 1 — y of the parent parton momentum pq. The frac- 
tional momentum y is defined with respect to the momenta carried by the colour connected 
particles: the quark, the antiquark and the gluon. In particular y is defined as, 

y = ^, (4.5) 

Vgqg 

and, since y qg > y m i n yqqg, the lower boundary of the y integral is y m i n . In this limit the four 
particle invariants are related to the invariants of the three remaining particles, 

Vqq = (1 - V) VQq, Vqg = V VQq, Vqqg = VQq ( 4 - 6 ) 

while the invariants containing p 7 become, 

y qi = {l-y) y Q -y, y ig = y y Ql . 

The matrix elements and phase space exhibit an overall factorization, 

l-H^rol 2 P q g^Q(y, s qg )\M Qqi \\ (4.7) 

with, |.M<2g 7 | 2 the three-particle matrix element squared for the scattering of a quark- 
antiquark pair with a photon, and, 



P q ^ Q (y,s qg ) = —P q ^ g (y). (4.8) 
s qg 

Similarly, the four particle phase space becomes, 

dPi d) (M;p q ,p g , Pl ,p g ) - dPi d \M;p Q ,p q ,p,)dPS(Pq,Pg,y) (4-9) 
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where dP 3 (M; pq, p q , p 7 ) is the three-particle Lorentz invariant phase space in d-dimensions. 
The collinear phase space factor reads ||15|| . 



(4tt) 



167r 2 r(l - e) 



(4.10) 



To evaluate the quark-gluon collinear factor, we need to integrate the collinear matrix 
element squared over the unresolved variables y qg and y. Reinserting the overall coupling 
factor, we have, 

&F{q) = J g 2 C F ^ 2e P qg ^ Q (y, s qg )dP^} (p q , p g , y) 

fa s C F \ ( 4nu 2 Y 2 ryminy q9g /-i 



1 /a s CV\ /47r/i 2 V 1 



2tt / ^ m 2 ; r(i -ef^yQ* 



2..- e (l-e)(4-6) T 2 (l- e 



-2/ 



e *mm 2e(l-2e) r(l - 2e) 

(4.11) 



We see that compared to the single quark-gluon collinear factor obtained considering the 
single quark-gluon collinear limit of a three parton cross section found in [15], the collinear 
factor given by eq. ( |4.11j) is multiplied by y^. This slight modification is caused by the 
change in the boundary of the y qg integration. The invariant y qg is here bounded to be less 
than y m iny q qg instead of y m \ n in the three parton case. 

Putting all the factors together, we find that in the collinear quark — gluon limit the four 
particle differential cross section da^ becomes, 

da 4 - C F (q) x (^j 2(2vr) 2 ^ J \M Qqi \ 2 dP^ (M;p Q ,p q , Pl ) = C F {q) x da Qqy , (4.12) 

where daQ qi is the three-particle differential cross section for the scattering of a quark- 
antiquark with an additional hard photon. All of the divergences are isolated in C F (q) and 
therefore d(jQ qi can be evaluated in 4-dimensions. 

In the region where the antiquark and the gluon are collinear the roles of quark and 
antiquark are exchanged. Now, the antiquark and gluon form a parent parton Q and the 
resulting contribution in this region of the four particle phase space therefore yields, 

dvqqyg -> C F {q) da qQ ^, (4.13) 

where, since y Qq = y q Q, 

C F {q) = C F {q). (4.14) 
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In order to match onto the single collinear quark-gluon regions, the soft gluon region is 
defined in eq. ( |3.7| ), 

Dqg *^ 2/min Dqqg i Dqg *^ 2/min Vqqg i Vq-f ^ 2/min j I/97 2/min • 

As in the collinear regions, the matrix elements and phase space both factorise in the soft 
gluon limit, 

\M mg 1 2 -> |A^| 2 / gf (<?), 

with the eikonal factor, 



(4.15) 



/«(*)- '' S "" 



SqgSqg 



and, 



dP 4 (d) (M;p 9 ,p 5 ,p 7 ,p 9 ) -> dP 3 (d) (M;p g ,p 5 ,p 7 )dP^ t (pg,p 5 ,p 9 ), 
where the soft phase space factor reads 111 



(4.16) 



(47T) 



167r 2 r(l - el 



'99 



S qg S qg 



'99 



As before, all of the dependence on the unresolved variables is collected into the soft approx- 
imations to the matrix elements and the phase space. We find, 



S, 



g 2 C F fi 2e f qq {g)dP { s ^ ft (p q , p q , p g ) 



a s Cp\ ( 4:7r fi 2 



2vr 



2 rym'myqq 



M 2 J r(l-e) y qq k 



9.9 



£/min£/qq 



99 



UqgUqg 



-e-1 



2tt 



1 2_ _ 2e _ e 

]Vf 2 i IY1 - e) e 2 VTa ^ y if 



(4.17) 



Again, apart from the changed boundaries, this is the same soft factor as given in eq. (3.33) 
of |l~5 |. As usual, the contribution to the cross section from the single soft singular region 
factorizes, 



da, 



9979 



Sp da qqi . 



(4.18) 



The sum of the single unresolved gluon contributions is then given by combining eqs. ( 4.12|) , 
(gig ) and ( gig) , 



da^ 79 -> (C F (q) + Cf(?) + da qqi = R^da, 



997' 



where the real unresolved factor R qqi depends on the invariant mass of the quark-antiquark 
pair and is given by, 



R, 



997 



a s C F 



2tt J Yil 



1 



' Ann 2 



x 



99 



H 2 ln 2 (j/ min ) - 3 \n(y qq y min ) + 7 



2tt 2 



(4.19) 
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This is the same real unresolved factor as defined in eq. (3.79) of |yj with y m - m — > y m - m yqq 
reflecting the changed boundaries of the unresolved gluon region of the four parton phase 
space compared with the unresolved gluon region of the three parton phase space. The 
singularities present in these single unresolved gluon contributions will cancel with those 
from the resolved photon one-loop 7* — > qqj process discussed in section |6T[ 



4.2.2 The unresolved photon contribution 



In the region where the quark and the photon are collinear we have, cf. eq. ( ft.4|) , 

Uqy ^ UrcYmUqq^ii Vq^i > VravaHqq^i V<ng - > 2/min; Vqg ^ 2/min- (4.20) 

The quark and the photon cluster to form a new parent parton Q such that each carries 
respectively a fraction z and 1 — z of the parent parton momentum pq. The four-particle 
matrix elements and phase space factorize in exactly the same way as in the quark-gluon 
collinear limit with the roles of photon and gluon being interchanged and y replaced by z. 
Unlike the quark-gluon case however, the photon is observed in the final state and hence 
only y qi is an unresolved variable. The momentum fraction z carried by the photon inside 
the quark-photon cluster is defined with respect to the momenta carried by the electromag- 
netically connected particles, 

z = ^. (4.21) 

Uqq"f 



In this limit, the four particle differential cross section again factorizes, 

da 4 -> C Fl (q)dz x daQqg (4.22) 

where daQ qg is the differential cross section for the production of a quark-antiquark pair and 
a gluon. The singular factor Cp 7 (q) is, 

c F M = - (^pj (jjjpj fTf^y 7 ^ VqI z ~ e (! - z )~ e p ^ z )- ( 4 - 23 ) 

There is a similar contribution where the photon is collinear with the antiquark. 



4.2.3 Overlapping of single collinear regions 

It is worth noting that in both, q — g and q — 7 collinear regions we have required y qig > y m \ n 
in order to guarantee that these regions match onto the double unresolved triple collinear 
region defined by y qig < y m i n . However this requirement has an important consequence. We 
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do not avoid the situation where both invariants y qi and y qg are simultaneously less than 
7/min and the constraints, 

Vq-/ ^ Vram Uqg *^ 2/min but Vqig ^* Umim (4.24) 

define an overlapping region of the two single collinear g — 7 and q — g regions. In this 
overlapping region, the matrix elements are correctly approximated by the sum of the two 
single collinear approximations and the error in the analytic contribution resulting from 
evaluating the approximated matrix elements over this restricted phase space region is of 
0(Vmm) an d therefore negligible. 

However, as we mentioned before, the photon + 1 jet rate is to be evaluated numerically 
applying the hybrid subtraction method. It is then crucial to ensure that the matrix elements 
are correctly approximated in this overlapping region since here we evaluate the difference 
between this approximation given by the sum of both collinear approximations and the 
full matrix elements. It is precisely to take into account these numerical contributions 
(which generate terms proportional to ln(y min )) that we are constrained to apply the hybrid 
subtraction scheme rather then the more commonly used phase space slicing method |jT4], ITH] 



where the two single collinear regions would have to be clearly distinct and not overlapping. 
However, in this region of phase space, either collinear limit on its own is a very poor 
approximation to the full matrix elements. As a consequence, using the phase space slicing 
approach, we would fail to obtain the necessary ln(7/ min ) cancellation in the physical photon 
+ 1 jet cross section. 



5 Double unresolved contributions 



In the previous section we have discussed the resolved and single unresolved real contributions 
relevant to the tree level process 7* — > qqjg. Each of these two classes of real contributions 
corresponds to final state configurations where more than two particles are theoretically 
"seen" and a 7 + 1 jet event can only arise if some final state particles are clustered together 
by the jet algorithm. Therefore, to allow the adjustment of our results to any jet algorithm 
used in the experimental analysis, the finite contributions to the differential cross section 
will be evaluated numerically. The two-particle unresolved contributions, on the other hand, 
already correspond to 7 + 1 jet events, and there is no further clustering by the jet algorithm 
needed. Because of this, the integrations can be performed analytically. 
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5.1 The triple collinear factor 



As we saw in section 3, the triple collinear contributions arise when the gluon and the photon 
are collinear to the quark. The triple collinear configuration is illustrated in Fig. ^.c.i. In 
order to evaluate these contributions we need to determine the appropriate approximations 
for the matrix element squared and phase space in the triple collinear limit and perform the 
phase space integrals over the unresolved variables. 

The triple collinear region of phase space is defined by eq. (|3.8f) , 

y<ri9 ^ Z/min and Vqg ^* 2/min- 

In this limit, y qgi is small and the photon, gluon and quark cluster to form a new parent 
parton Q such that, 

Pq+Pj+Pg =PQ- (5.1) 

The photon, the gluon and the quark carry respectively a fraction z, y and (1 — y — z) of 
the parent parton momentum pq, 

Pi = z p Q , P g = y Pq, p q = (l-z- y) p Q , (5.2) 
so that the invariants are given by the following, 

y q q = 0--y-z)vQq = 0--y-z), 

yq~< = zy Qq = z, (5.3) 

yqg = yygq = y- 

The algebraic structure of these double unresolved contributions is unique to the triple 
collinear limit of the matrix element squared, and when analytically integrated over the 
singular regions of phase space will form the triple collinear factor. These contributions 
are expected to arise in analytic calculations of exclusive quantities at the second order 
in perturbation theory. Such calculations have, to the best of our knowledge, not been 
performed before in the literature. 

We are interested in the triple collinear limit of the matrix element squared for the scat- 
tering of a quark-antiquark pair with a photon and a gluon. In this limit the d- dimensional 
four-particle matrix element squared factorises, 

l-Mgggcyl > Pggj^Q^Z, y, Sqry, Sqg, Sqgj) \-M-Qql , 

where lA^Qg] 2 is the two-particle matrix element squared and P qgi ^Q(z, y, s qi , s qg , s qgi ) de- 
fines the triple collinear splitting function. It is obtained by keeping only the terms containing 
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a pair of the unresolved invariants, s 97 , s qg and s qgi , in the denominator of the "full" four 
particle squared matrix elements \Ai qqgi \ 2 , 



+ 



4 (1 - z - y)(l + (1 - z - y) 2 - t{z 2 + zy + y 2 ) - e 2 zy) 

^qi^qg Z V 

4 {l-z-y)(l-z + e 2 zy) + (1 - yf - e(l - y){z 2 + zy + y 2 ) + e 2 zy 

^qi^qgi 

4 (1 - z - y)(l - y + e 2 zy) + (1 - z) 3 - e(l - ^)(z 2 + zy + y 2 ) + e 2 z?/ 
s qg s qgi Z V 
~ ^^f(l-6)^ + (l-e)^£-2eV (5.4) 

S qg~i \ S ig S 11 / 

This triple collinear splitting function is the generalization of the single collinear factor with 
three collinear particles instead of two and is as universal as the single soft and single collinear 
factors encountered earlier in section 4. 

In the triple collinear limit, we can reorganise the <i-dimensional four particle phase space 



given in eq. (4.2) into one part appropriate for the remnant Q — q pair multiplied by the 



integral over all of the unresolved variables, 

dP 4 (d) (M;p 9 ,p 7 ,p 9 ,Pg) -> dP^ d) (M;p Q ,Pq) x dP^ col (p Q ,p q ,p 1} p g ), 
where dP^ \M;pQ,p q ) is the usual two body phase space and, 

x dy qig dy q7 dy qg dy ig dz dy 8{y qi + y qg + y ig - y qig ). (5.5) 

Here the angular integrations for the rotation of the g^7-system around the q axis and for 
the parity of the g^7-system have been performed. In terms of the unresolved variables, the 
Gram determinant is, 

A 4 -»• ((1 - y - z) y ig - yy q7 - z y qg f - Azy y rj y qg . 

Reinserting all coupling factors into the matrix elements, we find that in the triple 
collinear limit, the cross section factorises, 



do"4 = 



a s C F \ ( ae 2 q 
2tt ) \ 2ir 



^4(27r)V e / \M qm \ 2 dPi d) (M;p q ,p q ,p„p g ) 



-> TC Fy dzx J\M Qq \ 2 dP^(M;p Q ,p q ) 

= TC Fy dz x a . (5.6) 
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As usual, (To is the two-particle cross section while the dimensionless factor TCp 7 dz contain- 
ing all the singularities is formally given by, 



TCF-ydz 



2vr 



-A, 



ae„ 



2tt 



'47T/i 2 



T(l - 2e) 4vr 



Pa. 



xdy qig dy qi dy qg dy ig dz dy 8(y qi + y qg + y ig - y qig ). 



(5.7) 



To evaluate TCp^dz, we must integrate out the unresolved variables over the phase space 
region where the triple collinear approximation is appropriate. This is, 



< y qig < y 



<y<l-z, 



while the other variables are constrained by the Gram determinant. Using the delta function 
to eliminate y ig , we find that the bounds on y qi are given by solving the quadratic equation 
A 4 = 0, which generates the additional constraint y qg < (1 — z)y qig . It is always most conve- 
nient to integrate the invariant mass that does not appear in P qgi ^Q first. This removes the 
Gram determinant and generates factors that regulate the other singularities. The first term 
in P qgi ^Q does not allow this approach and is rather more tricky, yielding Hypergeometric 
functions with complicated arguments after the first integration p7| . However, the integrals 
can be carried through and we find that after making a series expansion in e, the integrated 
triple collinear factor is given by, 
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6 3 y v ' \2 

+ \n(z) ln(l - z) (2 + ^ - 3 P^ 7 (z)) + ln(z)Li 2 (l - z) (4 - 2 z) 
+ ln 2 (l - z) (2 + ^ - 3 P g ( °! 7 (z)) + ln(l - z)Li 2 (l - z) (a - 2 z + 2 P<°1 



4 



- 7P(°) 7 (z) + ln 3 (z) ln 2 (^) ln(l - z) (3 - ^ + P^ 7 (z)) 

+ ln(z) ln 2 (l - z) (2 - z + 3 P^zj) + ln(l - z) (~ - ^ + 7 Pf 7 (z)) 

+Li 3 (l - z) (-4 + 2 z - 2 P^ 7 (z)) + S^l - 2) (2 - z - 3 P^z)) 

7 ln 3 (l - z)Pi%(z) , s , m , s , s 4 ln(l - zW 2 pW (z) 
+ , g 71 ^ + 4zLi 2 (l - z) + 9 P^ 7 (z)C(3) i 



+ \n(y min ) (Li 2 (l - z) (4 - 2 z + 2 P^%(z)) + ln(z) ln(l - z) (4 - 2 z + 2 P^; 

+ ln 2 (.) (3 - y " P i°^)) " 27r2P f 7(Z) + - *) (4 + z - 6 
+ K*) (2 + ^-3 P g ( °! 7 + 5 (ln(l - z)f P(%(z) + ~ - ^ + 7 P^fz)) 
+ ln 2 (i/ mm ) (ln(z) (2 - z - 5 P^z)) + 2 - ^ - 3 P^%(z) - ln(l - z)P^ 7 (z)) 



19 ln 3 (y min )p(°) 



" ■ (5. 



5.2 The soft /collinear factor 

The soft/collinear configuration arises when the photon is collinear to the quark and the 
gluon is soft as illustrated in Fig. [|.c.ii. As before, to isolate the singularities we need to 
determine the appropriate approximations for the matrix element squared and phase space in 
the soft/collinear limit and perform the phase space integrals over the unresolved variables. 



The soft / collinear region of phase space is defined by eq. ( |3.9| ) , 

and y qg < y 

min 

and yq g < y 

miri' 

In this limit, the photon and the quark cluster to form a new parent parton Q such that, 

P q + P 1 =P Q , (5.9) 
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while the energy of the gluon tends to (p g — > 0). The photon and the quark carry respec- 
tively a fraction z and (1 — z) of the parent parton momentum pg, 

p^ = zp Q , p q = (l-z)p Q , (5.10) 

and the invariants are given by the following, 

y qq = (1 - z) y Qq = (1 - z) and y qi = zy Qq = z. (5.11) 



Once again, the matrix elements factorise in the soft /collinear limit defined above, 

I-^^stI * ^qgi— >Q ' s 97> s i9i b qg-> b qgi) l-^Q<?l > 

As usual, | M.Qq\ 2 is the two-particle matrix element squared while P qq ^q l (z, s qi , s qg , s qg , s qgi ) 
defines the soft /collinear approximation to the squared matrix elements. This is obtained 
by setting y = in the triple collinear splitting function, P ?g7 _»Q, given in eq. ( |5.4| ) and is 
therefore, 

Pqgi-^Q ( z > S| 37> s qg-> "%> s qgi) = ~ " ~ ( (1 ~ z ) + ~ " ~ ) Pq-*ii z )- (5-12) 

b qi b qg b qg \ b qg-y J 



In the soft /collinear limit, we can again rewrite the phase space as the phase space for 
the Qq pair multiplied by an integral over the unresolved variables, 

dPjt d) (M;p q , Pj ,p g ,p q ) -> dP^ d) (M;p Q ,p q ) dP { s d) ft/col (p Q ,p q , Pl ,p g ), 

where, 

,r,f« / \ 1 1 M 6 /4t\ 2£ r A ._i_ e 

dP soft/col (PQ,Pq,P„P 9 ) = f^^^A^y{jp) [ " A4] 2 

x dy qjg dy qi dy qg dy ig dy qg dz 6(y qi + y qg + y ig - y qig )(5.13) 

This factor is similar to the triple collinear phase space factor given by eq. ( |5.1|) with y 
replaced by y 2 4- 

As y qig is unconstrained in this region of phase space, we choose to rewrite A 4 as a 
quadratic in y qig , and, when performing the phase space integrals, will first integrate over 
y qig . With the definitions of the invariants y qq and y qi in the soft /collinear region of phase 
space, 

~~ ^4 = (1 — Z) {yqygb ~ 2/<?79 ) (^TO ~ Vq^/ga), 

with y qyga ,b given by, 

yqfga,b = ~ I 2/<ry(l z ) Uqg Vq^Vqg i 2 ^Jy q ^y q ~ g 1 \j q ~ g Z I . 
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Once again, the cross section factorises, 

da 4 = {^f) ^4( 27 r)V e j \M qqgi \ 2 dPf {M ]Vq ,p q , Pl ,p g ) 
-> SC Fj dzx J\M Qq \ 2 dP^ d \M;p Q ,p q ) 

= SC Fj dz x cr - (5.14) 
The singular dimensionless factor SCp-ydz is given by, 



ft/col 



oc s C F \ ( ae 2 \ ( Att/j, 2 



2tt / V 2tt / \ M 2 J T(l - 2e) 4vr 
x y [-^V^ P qg j!l%\z,y qj ,y qg ,y qg ,y qgj ) 

xdy 979 dy 97 dy TO dy 79 y & dz 6(y qi + y qg + y ig - y qig ). (5.15) 

To evaluate the soft /collinear differential factor SCp-ydz we need to integrate P s °J^q \ given 
by eq. (|5.12|) over the soft / collinear phase space given by eq. ( |5.13D . 

The Gram determinant fixes the allowed range of y qig , while the other three unresolved 
variables are constrained to be less than y m { n . As expected, the y qig integral generates factors 
regulating the other phase space integrals and we find, 

*>- ~ m 1 ^ ( 4 i£T (^) (£) B >* ■ 



5.3 The double collinear factor 

The double single collinear region of phase space is defined by the following constraints on 
the invariants, 

Vq-y < 2/mim y qg < Z/miii) (5-17) 

with the additional requirement, 

y qg >ymin, (5.18) 

because the gluon is collinear to the antiquark but is not soft. This configuration occurs 
when the photon and the quark form a collinear pair simultaneously with the gluon and the 
antiquark being collinear and is illustrated in Fig. |2|.c.iii. 

In this limit, the photon and the quark cluster to form a new parent parton Q such that, 

P q + P 1 =P Q , (5.19) 
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while the gluon and the antiquark cluster to form a new parent parton, Q with, 

Pq+Pg=PQ, (5.20) 

The photon and the quark carry respectively a fraction z and (1 — z) of the parent parton 
momentum pq, 

Pi = z Pq, Pg = (1 - z) Pq- (5.21) 

whereas the gluon and the antiquark each carry a fraction y and 1—y of the parent momentum 
Pq such that, 

Pg = y PQ, Pg = (1 - V) PQ- (5-22) 
The invariants can be redefined as follows, 

y q q = (1 - y)(l - z) y QQ = (1 - y)(l - z) 

Vq 9 = vi l ~ z )yQQ = y( l - z ) 
2/7 = z ( 1 ~y)yQQ = *0--v) 

y J9 = yzy Q q = yz. (5.23) 



Using the redefinitions of the invariants given by eq. ([5.23 ), the four particle matrix 



element squared factorizes in the double single collinear limit as follows, 

\M q q Jg \ 2 — > P q7 _>Q.q g ^Q(z,y, S g7 , Sqg)\M Q Q\ 2 (5-24) 

with, 

Pq~i^Q;qg^Q~( Z -> Vi S 97' S qg) = Pq~1^Q( Z , S q-y)Pqg^Q{y, S qg)- (5.25) 

In other words, the double single collinear factor is the product of two single collinear factors. 

In this limit the four-particle phase space again factorises, 

dP 4 (d) (M; p q , p q , p 7 , p g ) -> dP 2 (d) (M; p , Pq) x dPji ;e (p 9 , p 5 , p 7 , p 9 ). (5.26) 

As before, the unresolved phase space factor contains integrals over five unresolved variables. 
Note that unlike in the soft /collinear region, y ig is precisely defined in eq. ( |5.23|) . Conse- 
quently the triple invariant y qig is also fixed, y qig = y qi + y qg + y ig = y. This appears to 
reduce the number of independent variables by one. However, a closer look enables us to 
assert that there is no inconsistency in this procedure. In fact, the boundaries of the y qig 
integration generated by the Gram determinant turn out to be y q -yg a ,b — y ± C(z/mm). Hence 
by replacing y qig by y in order to obtain the double single collinear matrix elements we make 
an error of 0(y m i n ), which we do throughout the calculation and it is therefore a consistent 
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approximation to make. Integrating out y qig and the unresolved angular variables, we find 
that, 



,r,fcn / \ 1 /4vr\ 2e M 4 . . 
dP doLbie(P^ Pg, Pi, Pa) = ZT TH TTu^ITa d Vw dz z(l- z)y 



T 2 (l -e) \M 2 J 16(2tt) 4 
x(dy qg dy\y(l - y)y qg ] j , (5.27) 



which is exactly the product of two single collinear phase space factors as one could have 
expected. 

As with the previous double unresolved contributions, the integration of the resolved 
two particle matrix elements over the two particle phase space yields a factor of cr . This is 
multiplied by the integral of the approximation of the unresolved matrix elements over the 
unresolved phase space. Explicitly, we have, 

da qqig — > DC Fy dz x <t , (5.28) 

where, 

DC Fl dz = (^) (^) 4(2vr) 4 {^f J dP^ ble P q ^ Q . qg ^ Q (z, y, s q „ s qg ). (5.29) 

The unresolved region is specified by eq. (|3.10|) and the constraint y qg > y m i n fixes the lower 
boundary of the y integral to be y TO i n /(l — z) since in the double single collinear region 
Vqg = 2/(1 — z )- The integrals are straightforward, and we find, 



^\ 2e / q B C F \ ( ae 2 q 
l (l-e) 2 V M 2 J V 2tt / V2tt / 

x ^ f 2 ^ (1 _ _ (1 -'')(*- ') I^lziA . (5 .30) 

e 3 I yminV ; 2(1 - 2e) T(l - 2e) / V ; 



5.4 Strong Ordering 

As a check of our calculation of the real two particle unresolved contributions to the differ- 
ential cross section, we have rederived them in the strongly ordered limits | 27| . Instead of 
considering particle 1 and particle 2 to be collinear at the same time to particle 3, we con- 
sider the two different contributions; either particle 1 is collinear to particle 3 followed by 
particle 2 collinear to the cluster of particles 1 and 3, (denoted by (13)), so that (2/13 <C 7/23), 
or particle 2 is collinear to particle 3 followed by particle 1 being collinear to particle (23) 
where we have (2/23 *C 2/13). In general, within the strongly ordered approximation, each of 
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the unresolved real contributions (triple collinear, soft/collinear and double single collinear), 
gets "replaced" by the sum of two strongly ordered contributions. In addition to changing 
the approximations to the matrix elements, the phase space is also reorganized to be the 
product of two single unresolved factors. However, although the strongly ordered approxima- 
tion correctly reproduces the leading divergences -those proportional to 0(1/ e 3 ) or 0(1/ e 2 ) 
which are associated with the leading and next-to-leading logarithms - it does not generate 
the correct subleading divergences proportional to 0(1/ e). The finite terms of 0(1) are also 
incorrectly reproduced. We understand this as follows: The poles in 1/e 2 , 1/e 3 arise from 
the evaluation of successive phase space integrals at the lower boundaries where the strongly 
ordered approximation is very close to the "full" approximation of the matrix elements. On 
the other hand, terms proportional to 1/e arise when evaluating only one of the phase space 
integrals at its lower boundary while the other phase space integrals contain significant con- 
tributions close to their upper boundaries. At these upper boundaries, the two invariants 
defining the strongly ordered limit are no longer strongly ordered and the approximation is 
invalid. 



6 Virtual contributions 



In the previous two sections we have decomposed the four-particle phase space and ex- 
tracted the divergences present in the 0(aa s ) four-parton process 7* — > qcpjg where one or 
two particles are theoretically unresolved. In other words, only two or three particles are 
theoretically identified in the final state. If three particles are theoretically well separated, 
the experimental cuts will combine these particles further to select photon + 1 jet events. 
In this section we will take into account the exchange of a virtual gluon in the 7* — > qq^ 
process, which when interfered with the tree level process also gives rise to contributions of 
0(aa s ). 

As discussed in section [|, the calculation naturally divides into two parts, depending on 
whether or not the three particles are resolved. Both resolved and unresolved contributions 
are divergent and need to be combined with the appropriate real contributions described 
earlier in sections 4 and 5. For the resolved virtual contribution, the quark, the antiquark 
and the photon are clearly distinguishable and the divergences will cancel when combined 
with the real contributions where the gluon is either collinear with one of the quarks or is soft 
and the photon is theoretically resolved (c.f. section [|). On the other hand, in the unresolved 
virtual contributions, the quark and photon are collinear and form a single pseudo particle, 
Q, the parent quark. The leading singularity occurs from a soft gluon being internally 
exchanged simultaneously with the collinear emission of the photon from a quark and is 
proportional to P q ^(z)/e 3 . These most singular poles will cancel with those present in the 
soft/collinear contributions calculated in the previous section. 
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6.1 The resolved contribution 



The "squared" matrix elements arising from the 7* — > qq'f process at one loop is part of the 
0(a 2 ) corrections to the three-jet rate in e + e~ annihilation, which was originally derived 
by Ellis, Ross and Terrano in P2| . As we are interested in the virtual contributions with 
an outgoing photon instead of an outgoing gluon, we need to replace the colour factors in 
eq. (2.20) of @ with, 

C A -> 0, Np — > 0, C 2 F -> C F , 



as well as the replacement, 
when the quark has charge e q . The contribution to the cross section can be written as, 



a 2 — >■ a s ae 2 , 
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(6.1) 
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The function R is defined as, 



R(yqq,yq~,) + 



2 In y qq 



Vqq + (2/99- + Viftf 
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iqq 
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99 



(y qi + y^) 2 Vq'y t Vq-y ) J 



.. (6.4) 



R(x,y) = lnxlny — lnxln(l — x) — lnylnfl — y) H — tt 2 — Li 2 (x) — Li 2 (y). (6.5) 

6 
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To ensure that the photon is resolved from the quark and antiquark, we define the resolved 
three parton phase space to be, 

Vq-y ^ 2/miri) Dqi ^ 2/min- 

In this region the resolved virtual cross section can be written as, 

d(J qqi = V 1^ d(T Sll + da f?7' ( 6 - 6 ) 

with, 

(cte^\ 1 r 

2^ J 2 ^ 2 J^p J F (fe%7» Vqi) dP 3(M;p g ,Pg,Py). (6.7) 

The singularities present in V q q 1 precisely cancel with those from the 7* — > qq'jg process when 
the gluon is unresolved given in eq. ( [4.19| ). The contribution from the finite function do"^ 7 
will be evaluated numerically using the expression of Ellis, Ross and Terrano in eq. (|6.4|) 
and the experimental jet algorithm to select photon + 1 jet final state events. The lowest 
order resolved cross section dcr^ 7 will also be evaluated numerically. The more interesting 
problem lies in the unresolved region as we shall see in the next subsection. 



6.2 The unresolved contribution 

In the unresolved region of phase space, the quark becomes collinear with the photon so that, 
defining a new parent quark, Q, with momentum pq we have, 

Pq+P~, =PQ- 

As usual, we introduce the variable z, 

p q = (l-z)p Q , p 1 = zp Q . (6.8) 

The photon carries then a fraction z of the composite quark momentum. In this single 
collinear limit the three particle phase space factorizes into a single collinear phase space 
factor, as we saw in section |], 

dPf(M; M ,p 7 ) - dP^{M; PQ ,p q )dP^(p q , Pl ,z). (6.9) 

At this stage we would like to take the corresponding limit of the virtual matrix elements. 
However, we note that the form given in eq. (|6.4| ) is unsuitable for taking the collinear limit, 
since as y qi — > 0, terms of the form, 

log {Vgi) 
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are present. Such terms are problematic and are generated by taking the y qy — > limit after 
an expansion of the virtual matrix elements as a series in e. The correct procedure is to take 
the collinear limit first and then expand the matrix elements as power series in e. Doing this, 
we find that the matrix elements factorize and are proportional to the lowest order 7* — > Qq 
matrix elements, 

^V col (z,s qi )\M Qq \ 2 , (6.10) 

where, 

v ^ s ^ ~ [ms) r(l-2e) 

( ( 2 2 2 1 (3 + 2e) 

x \P q _ n (z)[ — - + — y~ e - —y~*(l-z)- £ F 21 {-e,-e;l-e;z^ 



e 2 1 e 2 tf 97 e 2 y 97 ^ > ZL V ' ' ' "/ e (1 — 2e) 

ez - 1) 
1 -2e) 



The Hypergeometric function F 2 i(a, b; c; z) can be expanded as a series in e, 
F 21 (-e, -e; 1 - e; z) = 1 + e 2 Li 2 (z) + e 3 (Li 3 (z) - S 12 (z)) + G(e 4 ) 



The collinear behaviour of one-loop amplitudes has been studied by Bern, Dixon, Dunbar 
and Kosower |33| using helicity amplitudes. We have checked our result for the collinear limit 
of the squared matrix element by comparing with their results. 

Integrating the unresolved squared matrix element over the corresponding unresolved 
phase space region yields a dimensionless virtual collinear factor VCpydz multiplying the 
lowest order two particle cross section, 



dc^ 7 -> VC Fl dz x (T , (6.12) 



where, 



/in m 



Performing the y qi integration, we find, 



(6.13) 



27r J \ 2tt J \ M 2 J T(l - 2e 
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+ <;-">). (6.i4) 
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As expected, the most divergent term in this expression is proportional to P q ^(z) and 
precisely cancels the leading singularity present in the two-particle unresolved contributions 
from the four parton process discussed in the previous section, namely the leading singular- 
ity in the soft /collinear contribution SCf^ (c.f. section |5.2| ). The subleading poles do not 
cancel and are ultimately factorized into the 0(aa s ) counterterm of the quark-to-photon 
fragmentation function, as will be discussed in section ||. 



7 Fragmentation contributions 

In addition to the real and virtual contributions derived in the three preceding sections, we 
need to consider contributions associated with the fragmentation processes shown in Figs. |l|.c 
and |l|.d. These arise when a quark-antiquark pair associated with a real or virtual gluon are 
produced, followed by the fragmentation of a quark into a photon. The contribution of these 
processes to the differential photon + 1 jet cross section is determined by the convolution of 
the tree level 7* — > qqg or the one-loop 7* — > qq cross section with the bare quark-to-photon 
fragmentation function, (x), and has the symbolic form, 

da %j(g) = dcT OT(9) D f-^y( x ) dx - ( 7A ) 

Here, dcr^ s ), and do- gq -( g ) are the fully differential cross sections and x is the ratio between 
the photon and the parent quark momenta. 

The bare quark-to-photon fragmentation function D?^{x) is the sum of a non per- 
turbative part, D g ^ 7 (x, /zp) which depends on the factorization scale /zp and can only be 
determined by experiment, and a perturbative counterterm. Since the underlying 7* — > qq(g) 
process is already of 0(a s ), only the 0(a) counterterm needs to be considered and we have 
cf. eq. (F21D , 

Bf-» = IW*. „,) + \ (|£) )' j^L_i$»,(x). (7.2) 

As usual, this separation introduces a dependence on the fragmentation scale /Zp in both the 
counterterm and the physical fragmentation function D q ^(x, /zp). 

As discussed in section [3], the fragmentation contributions separate into three categories, 
depending whether the gluon is resolved, unresolved or virtual. If the gluon is identified 
in the final state, we will find that the singularities present in this resolved fragmentation 
contribution are exactly cancelled by the real collinear photon/resolved gluon contribution 
from the 7* — > qqg'j process. On the other hand, if the gluon is unresolved, it may be 
collinear with the quark or with the antiquark or it may be soft. In the absence of the 
quark-to-photon fragmentation function, the infrared singularities from the 7* — > qqg process 



37 



exactly cancel against those from the one-loop 7* — > qq process. However, because of the 
fragmentation function, this is no longer the case. When the gluon is collinear to the quark 
which subsequently fragments into a photon, the parent quark momentum is shared between 
the quark and the gluon and the fractional momenta carried by the photon and the gluon 
are related to each other. This introduces a convolution between fragmentation function 
and parton level cross section. As we shall see in section |7T^ , a large part of the divergences 
present in this contribution cancels against the singularities present in the double unresolved 
contributions discussed in section |[ 



7.1 Resolved contributions 

To ensure that the gluon is resolved from the quark and antiquark, we define the resolved 
three parton phase space to be, 

Vqg ^ l/min; Uqg Umm- 

Events satisfying this constraint contribute to the 7 + 1 jet differential cross section in the 
following two cases: 

(i) The gluon is clustered together with the quark which fragments into a photon. 

(ii) The gluon is clustered to the antiquark or it is isolated while the antiquark is clustered 
with the photon jet. 

In both cases the cross section has the form given by eq. ( [7. 1| ) with x, the fractional mo- 
mentum carried by the photon inside the quark-photon collinear cluster. Note that x is a 
theoretical parameter which is only related to the momenta of quark and photon. It does not 
necessarily coincide with the fractional momenta carried by the photon inside the photon 
jet, z, which is reconstructed by the jet algorithm. In particular x = z only holds if the 
photon jet only contains the quark and photon, while the antiquark and gluon are combined 
to form the second jet. If on the other hand, the antiquark or the gluon are clustered by 
the jet algorithm into the photon jet, one will generally find z < x. Ultimately, it is the 
experimental z, which is required to be greater than the experimental cut z cut . 

We note that the singularity structure from the qqg^f final state in the limit where the 
quark and photon are collinear (discussed in section |4 . 2 . 2|) is proportional to P q ^(x) and 



depends only on the theoretical x value. In fact, when the gluon is resolved, the cancellation 
of the singularities between the qqg final state with fragmentation counterterm and those 
generated in the qqgj final state when the quark and photon are collinear is unaffected by 
the possible discrepancy between x and z. This explicit cancellation will be demonstrated 
in section R. 
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7.2 Unresolved contributions 



In the previous subsection, the precise value of z was determined by the jet algorithm and 
is not necessarily the same as x. Similarly, when the gluon is unresolved, z and x do not 
necessarily coincide. If the gluon is virtual, soft or collinear to the antiquark, we can identify 
the ratio between the photon and the quark momenta x by z, since only fragmenting quark 
and photon form the "photon" jet. and the cross section has the form given by eq. (|7. 1|) . In 
each case the partonic cross section da qqg factorizes into a single unresolved factor multiplying 
the tree level cross section cr - As these factors have already been derived before flS] , we 
will merely quote their results. If a gluon is exchanged internally, the contribution to the 
7* — > 7 + 1 jet rate reads, 

da° = a V qg Df^{z)dz, (7.3) 

with, 

a s C F \ /47r/i 2 yr(l + e)r 2 (l-e) , 2 3 Q , 2 ia , 3 2 2 , 



v - = [-if) [-w ) r ( i- 2e ) [-7 - s - 8 + ^ - 166 + r e + 

(7.4) 

When the gluon is real but soft, the invariants y qg and y qg are both less than the theoretical 
cut y min we find, 

daf = a S F Df^(z)dz, (7.5) 

with the soft gluon factor being, 



2tt J \ M 2 J T(l - e) V e 2 
Similarly, when the gluon is collinear to the antiquark, y qg < y min but y qg > y m i n , we have, 

d<,) = 0(9)^(^(12, (7.7) 

with, 

r (^-( asCF \( A7[ ^\ 1 f 2 | (l-e)(4-e) T 2 (l-6) \ 

™ V 2tt / I M 2 J r(l-e) ymin l e 2 ymin 2e 2 (l - 2e) T(l - 2e) /' K ' 



On the other hand, if the gluon is collinear to the quark, so that the gluon carries a 
fraction y of the quark/gluon cluster momentum, z is no longer equal to x. In fact, z is 
given by the product of the momentum fraction carried by the quark, 1 — y and the ratio of 
the photon and quark momenta x, so that, 

z = x(l - y). 
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We therefore introduce the constraint, Jq dz 5 (x(l — y) — z) and integrate over x so that, in 
the collinear region, 

V qg < 2/min and y qg = y > y min , 

we have, 



da c(g) 



-a 



x 



dy 



l-z 

Vmin l — y 



v M 2 J T(l - e) 
[y(l-y)rP q ^ g (y)D 



B 

<?^7 



dz. 



(7.9) 



The trivial integral over y qg has been carried out while the factor of (1 — y) comes from the x 
integration. The y integral now involves the fragmentation function and requires some work 
to evaluate. The resulting expression will involve a convolution of the splitting function with 
the fragmentation function. However, the convolution integral present in eq. (|7.9|) appears to 
have an explicit y m [ n dependence coming from the lower boundary of the y integral. We know 
that since y m ; n is an artificial parameter which cannot influence the physical cross section 
for any choice of fragmentation function, the y m i n dependence must only act multiplicatively 
on the fragmentation function. 

We therefore add and subtract the contribution where a gluon is collinear to a quark 
multiplied by Dg_ n (z). The convolution integral can be rewritten in the following way, 



dtT C(q) + 



a s C F 
2tt 



dy [y(l - yT e P q ^ g (y) 



r(l-e) 
2 , 



Vmtn Df_^{z)dz 

_ (l-e)(4-e) T 2 (l-e) 
2e(l-2e) T(l - 2e) 



a C F (q)D"Jz)dz 



so that C F (q) = C F (q) and with da®, y given by, 



(7.10) 
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dy[y(l-y)]- e P g ^ g {y) 
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1-2/ 



dy [ y (l- y )]^P q ^ g (y)D 



q^i v 



In the first integral in the expression for da^ q y, the integrand vanishes when y — > 0, and we 
can safely extend the range of integration to 0. By doing so, the convolution contribution 
itself becomes y m i n independent as required. 
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Recalling the definition of the Altarelli-Parisi splitting function |24 

1 + (1 - yf - ey 2 



Wi/) 



y 



making the change of variables y = 1 — t, and using the "+" -prescription to evaluate the 
singular parts of the integrals, we can recast the expression for da^ q y in a more suggestive 



form, 



da E(iY 
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(7.11) 



As expected, the coefficient of the leading pole term is the universal lowest order Altarelli- 
Parisi |23| quark-to-quark splitting function in four dimensions, see eq. fl2.17| ). 



The contributions directly proportional to given in eqs. ([7.4D , (|7.6| ), (|7.8|) and by 



the second term of eq. ( |7.10| ) can be combined together and yield a finite result, 



da D K 



a (C F (q) + C F (q) + S F + V qq ) Df^(z)dz 



where K, qq is the finite two quark /C-factor given in eq. (4.22) of Ref. |l5fl . Since the 0(a) 
bare fragmentation function counterterm is proportional to 1/e, we need to expand K qq up 

to 0(6), 
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(7.12) 
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Finally the sum of all unresolved contributions involving can be written in terms 

of the dimensionless fragmentation collinear factor FC Fj , 

d<T qq(g) = da K + da c( q y 

= FC Fl dz x cr . (7.13) 

Expanding up to 0(e), we find, 

/ A„.,2\ e ^ A+ /.vN r 1 

.(2) 



- rn^T (^) (#) ' I', T ^ Ci) - [-7^« + ^ + 



(7.14) 



where, 
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2 \ l-t ) + 21- t y 1 l-t y ' 

+ (1 -t)ln(t) + (1 -t)ln(l -t). (7.16) 



The final step is to insert the decomposition of the bare fragmentation function given 
in eq. ( |7.2| ) into the sum of all fragmentation contributions given in eq. ( |7. 14j ) . We divide 
the fragmentation collinear factor into two pieces, one depending on the non-perturbative 
fragmentation function and one involving the perturbative counterterm, 



(7.17) 



with, 



FC p Fl 




D 



9^7 



Z, fJ, F ), 



-Lp(0) +c (l) +ec (2) 



Ip(0) 

(7.18) 



We see that this unresolved fragmentation contribution contains 1/e 2 poles as leading sin- 
gularities which must be combined with the virtual and the double unresolved singularities 
presented in sections || and |5|. 
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8 Sum of all contributions 



So far we have isolated all the divergent contributions to the 7 + 1 jet rate arising in the 
unresolved regions of the phase space analytically. However, to obtain a physical (finite) 
next-to-leading order prediction for the photon + 1 jet rate, we must combine these di- 
vergent contributions together and absorb the remaining singularities into the perturbative 
counterterm of the bare quark-to-photon fragmentation function. Following a fixed order 
approach to evaluate the photon + 1 jet differential cross section, we evaluate this rate order 
by order and consequently determine the perturbative counterterm of the bare fragmenta- 
tion function order by order. The 0(a) counterterm has been determined while performing 
the calculation up to order a |Q, we shall now utilise the 0(aa s ) perturbative counterterm 



given in eq. (|2.21|). 



Recalling the generic structure for the photon + 1 jet rate given in eq. ( |2.1| ) we have, 

+ -^®D* (8.1) 



1 dgfljet+V) _ 1 da, | 2 da g ^ DB 
<tq dz <r dz do dz y ~ 



Here, the factor of 2 reflects the fact that by charge conjugation, the quark and antiquark 
fragmentation contributions are equal. We recall that at 0(a), the differential cross section 
is given by, 

1 da LO (l iet + V) 1 da R - 

— 1 J = 2X> + — — (8.2) 

cr dz (T dz 

where V q ^ is the effective fragmentation function of eq. ( j2.13p [f|. Then at 0(aa s ), da 7 
receives contributions from the tree level four parton process and the one-loop three parton 
process. In sections [| and §, the singularities in the four parton process were isolated so 
that, 

1 f\rr - 1 diT^ 1 dfT^ drr^ 

+ iX g g 7 h Z O^-y + OOf 7 + ULyp^j + ZOf 7 . ^O. 6) 



(T dz (T dz (T dz cr 

The divergences are concentrated in factors multiplying finite resolved parton differential 
cross sections. Similarly, the singularities from the one-loop process were reorganized in 
section ^|, 

= V„.,^ + ^ + 2^ FT . (8.4) 
(T dz (To dz (T dz 

The fragmentation contribution from the real 7* — > qqg and 7* — > gg processes discussed in 
section [?| can be written, 

1 daq ' (9) ~ D* = -do*D*{z) + 2FC FT (8.5) 



(T dz 9->7 (J 
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Finally, there is a contribution from the 0(aa s ) perturbative counterterm in -D^L 7 multi- 
plying the lowest order 7* — > qq cross section, cf. eq. ( |2.21|) , 

v^® D ^ = 2D ^ {z) - (8 - 6) 



The overall result for the photon + 1 jet cross section is obtained by summing these four 
contributions. Regrouping the 0(aa s ) terms according to the resolved parton cross section 
they are proportional to, we obtain, 



1 da NLO (l jet + V 
cr dz 



1 da^ 1 dai 
\z 
da? 



_i_ir _ Z _332 _i_ 111 
cr dz a dz 



+2Vq ^ — 22£ 
^0 

1 da R - 

+±^ni. ( 8 . 7 ) 
cr d2 

Recalling the definition of C Fl together with the 0(a) perturbative counterterm given in 
eq. (|2.13|) yields the effective quark-to-photon fragmentation function up to O(a), 



B 



while, 



,^ 2 log (y mm )-31og( 2/mm ^-) + T -lj. (8.9) 

These, and the other contributions given in the last three lines of eq. (|8.7| ) are finite and 
will be evaluated numerically. The sum of the double unresolved factors, together with the 
fragmentation counterterm is also finite and provides the 0(aa s ) contribution to the effective 
fragmentation function, 
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where is given by eq. ( |7.16| ). 

Note that although each contribution given in eq. ( |8.7| ) is now finite the individual con- 
tributions depend on both the factorization scale /xp and the parton resolution parameter 
2/min- When combined numerically this y m \ n dependence must cancel as we will explicitly 
show in the next section. 



8.1 The evolution equation for D 



/i F ) 



In order to obtain a finite differential cross section we have factorized the collinear singu- 
larities in the perturbative counterterm of the bare quark-to-photon fragmentation function 
at some factorization scale /xp. The bare quark-to-photon fragmentation function should 
however not depend on the scale at which the factorization procedure takes place. Requir- 
ing that it is independent of the factorization scale /xp yields the evolution equation for 
the renormalized non-perturbative fragmentation function D q ^(z, /ip). We can see this by 
differentiating D q with respect to ln(/xfi), 



and using eq. ( |2.21|) , 
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For terms in the third line of this equation, which are proportional to a s , the variation of 
the non-perturbative fragmentation function with respect to /xp, dD q _+ 7 (z, /Xp)/S In (/Xp) is 
given by the lowest order evolution equation for D q _,^(z, /ip). To be more precise, at 0(a) 
and keeping terms of 0(e), we have, 
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Inserting this into eq. flS.12|) , we find the expected, \i independent, result, 
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8.2 An exact solution of the evolution equation 



The evolution equation ( 8.14j ) is insufficient to uniquely determine the non-perturbative 
quark-to-photon fragmentation function D q ^(z, /ip). However it is possible to give an exact 
(up to 0(aa s )) solution of this next-to- leading order evolution equation. This solution is a 
first step leading to the ultimate determination of Dg^lz, fip)- In the same way, the exact 
(up to 0(a)) solution of the leading order evolution equation [f| led to a determination of 
the quark-to-photon fragmentation function by a comparison between the LO calculation of 
the photon + 1 jet rate and the data [§]. 

We construct the exact 0(aa s ) solution by imposing that it takes the following general 
form, 
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where A, B are unknown functions of z, and the constant of integration //q. Here, 
D^~J(z, ix f) is the exact solution of the lowest order evolution equation obtained by ignoring 
the terms proportional to a s in eq. ( fg,14|) , 



(8.16) 



where the non-perturbative input fragmentation function D q ^(z,fi ) is to be determined 
by data. Inserting eq. ( |8.16| ) in the general form for the exact solution of the next-to-leading 
evolution equation and neglecting all terms which have more than one power of a s , we obtain, 



A = Pi ! In 



q^q ■ 



B 



P (± : ln(^)- 1 -P(% ® Pj°L In 2 ( & 



9^7 ' 



9^7 ■ 



Mo 



(8.17) 



so that the exact solution of the next-to-leading order evolution equation reads, 
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This solution has some interesting properties. First, it is exact at the order of the 
calculation (0(aa s )), and yields no terms of higher orders. The photon + 1 jet rate with 
this solution implemented is therefore independent of the choice of the factorization scale 
fiF, as it should be in a calculation at fixed order in perturbation theory. Second, we 
have not resummed terms proportional to In /if*. We choose to do this for a variety of 
reasons. Such resummations are only unambigous if the resummed logarithm is the only 
large logarithm in the kinematical region under consideration. If logarithms of different 
arguments can become simultaneously large, resummation of one of these logarithms at a 
given order implies that all other potentially large logarithms are shifted into a higher order 
of the perturbative expansion, i.e. have to be neglected. In our calculation, we encounter 
three different potentially large logarithms, In up, ln(l — z) and ln?/ cu t, and resummation of 
one of these would immediately imply that effects from the other logarithms at the given 
order of the calculation are neglected, a procedure which appears to be clearly inconsistent. 
Furthermore, since the overall photon + 1 jet rate does not depend on /ip, it appears to 
be more than doubtful that such logarithms should be resummed. As a final point, while 
Hf is an artificial parameter, it does represent the boundary between the perturbative and 
non-perturbative contributions. It corresponds approximately to the transverse momentum 
of the photon with respect to the jet axis, a physical scale, typically of the order of a few 
GeV, and not a large scale such as the e + e~ centre-of-mass energy. 



9 Dependence on y 



mm 



We have now collected all necessary ingredients to evaluate the 0(aa s ) 7 + 1 jet differen- 
tial cross section numerically. There are four separate contributions specified by eq. (|8.7|) 
and determined according to the number of partons in the final state and the presence 
or absence of the quark-to-photon fragmentation function. For each contribution, the ap- 
propriate matrix elements are integrated over the corresponding phase space using Monte 



Carlo techniques, i.e. with VEGAS [34]. In particular, for each event the invariants yij are 
defined allowing the reconstruction of the four-momenta of the particles in the events. 
The Durham jet clustering algorithm, with a particular jet resolution parameter y cu t, is then 
applied to these momenta to select 7 + 1 jet events. One of the clusters is identified as a 
photon if the fraction of electromagnetic energy z inside the jet is greater than z cut , a value 
fixed experimentally. Moreover, the fragmentation function mentioned above is the effective 



fragmentation function which at 0(a) is given by eq. ( [2.21|) and at 0(aa s ) is the sum of 



the fip dependent quark-to-photon fragmentation function and a finite contribution coming 



from the unresolved photon contributions given by eq. ( |8.10| ). The four contributions are: 



(A) 2 partons + photon 

This contribution corresponds to the process 7* — > qq'j with a hard photon in the 
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final state. It is present at 0(a) and 0(aa s ) (due to the presence of a theoretically 
unresolved real or virtual gluon). 

(B) 2 partons + "fragmentation" 

The 7* — > qqf&Dq^ process is present at 0(a) and 0(aa s ). In particular, it contains 
the finite terms corresponding to the q — 7 collinear region at 0(a) and the finite parts 
associated with all double unresolved regions at 0(aa s ). 

(C) 3 partons + photon 

This contribution is only present at 0(aa s ) and describes the 7* — > qqjg process where 
both photon and gluon are theoretically resolved. 

(D) 3 partons + "fragmentation" 

The 7* — > qqg ® D q ^ process with a hard gluon in the final state, is only present at 
0(aa s ). 

Each of the contributions depends on the theoretical parton resolution parameter 2/ mm . 
In Fig. I we show the various contributions to the integrated cross section 0.7 < z < 1, 
for a single quark of unit charge. The jet resolution parameter is y cut = 0.1 and only the 
0(aa s ) contributions are shown. As can be seen in Fig. |, the magnitude of the various 
contributions to the differential cross section increases dramatically as y m \ n becomes smaller. 
This rapid rise is due to the presence of logarithms of y m i n due to expanding the residues of 
the poles in e. Since the leading poles are (9(l/e 3 ), powers of logarithms up to ln 3 (y mm ) are 
present. 

The final integrated cross section, which is the sum of all theoretically resolved and 
unresolved contributions, must of course be independent of the artificial parameter 2/ mm . Each 
individual term has a very strong dependence on y m i n , but, as can be seen in Fig. ||, the sum 
of all resolved and unresolved contributions is clearly y min independent (within the numerical 
errors of the calculation) provided y m i n is taken small enough. In practice, this means values 
of ymin ranging between 10~ 5 to 10~ 9 for the chosen value of the experimental jet resolution 
parameter y cnt = 0.1. For y min = 10~ 8 , the magnitude of the individual terms is (9(5000), 
while the final result (after enormous cancellations) is O(10). This figure demonstrates the 
consistency of our approach - there is a region of parameter space where the choice of the 
unphysical parameter y min does not affect the physically observable cross section. Note that 
for large values of y min the cross section deviates from the y min -independent value. This is 
because for large y min values the approximations used in the analytic calculation become 
less accurate. In particular, terms of 0(y m - m ln 2 (y min )), which have been neglected, become 
sizeable. On the other hand, for values of y m i n below 10~ 9 the errors on the result become 
important due to the necessity of cancelling large logarithms numerically. The overall result 
becomes therefore less stable numerically for such small values of y mm . A reasonable choice of 
?/ m i n , which does not lead to problems of numerical instability is therefore y m i n /y C ut = 10~ 5 . 
We shall use this value when comparing our results with the ALEPH data in the next section. 
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Figure 5: Contributions of the individual terms ((A) , (B) , (C) , (D)) to the total cross section 
as function of y min for y cut = 0.1 and z cut = 0.7. For clarity, only the next-to-leading order 
contributions are shown. Furthermore we take ae 2 q = 2tt and a s Cp = 27r. 



10 Results 



In the previous section, we have demonstrated the consistency of our approach to evaluate 
the photon + 1 jet rate at 0(aa s ) by showing that the results were y min independent. It 
is therefore possible now to use these results to determine the non-perturbative quark-to- 
photon fragmentation function D q ^(z, fip) up to this order from a comparison with the 
experimental data || from the ALEPH Collaboration at CERN. As a first application we 
will use this newly determined fragmentation function to evaluate the isolated photon + 1 
jet rate in the democratic clustering approach as a function of the jet resolution parameter 
y cnt . These will also be compared to experimental data ||. 



50 



— 14 

N 




Figure 6: The sum of all next-to-leading order contributions to the total cross section as 
function of y min for y cut = 0.1 and z cut = 0.7. Only the next-to-leading order corrections 
are shown and we take ae 2 q = 2tt and a s Cp = 27r. The solid line is a fit of the form 

a + fr/min In 2 2/min- 



10.1 A determination of D q ^(z, fip) at next-to-leading order 



The quark-to-photon fragmentation function D q ^(z, fip) which appears explicitly in the 
expression of the photon + 1 jet rate given in eq. ( |8.7| ) has been given as the exact solution of 
the perturbative 0(aa s ) evolution equation by eq. (|8.18|) . As mentioned before, all unknown 
non-perturbative contributions to this function are contained in its initial value, D™p(z, /io) 
which needs to be extracted from the data. We perform a three parameter fit (/io is fitted as 
well) to the ALEPH 'photon' + 1 jet data |J for the z distribution, at a jet resolution 



parameter y cut = 0.06 and a s (M z ) = 0.124 and obtain [[12 




D^ LO \z,fx )= -r-M (~if] 7 (z)\n(l-z) 2 + 20.8(1 -*) - 11.07), (10.1 
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where /io = 0.64 GeV. For reference, the lowest order fit obtained by the ALEPH Collabo- 
ration |J is, 

DZ%°\z,lM>)= (^) (-P^)ln(l-,) 2 - 13.26), (10.2) 

with /io = 0.14 GeV. In both cases, the logarithmic term proportional to P^lJz) is intro- 
duced to ensure that the predicted z distribution is well behaved as z — ► 1 Q. 

The results of the 0(aa s ) calculation using this fitted next-to-leading order fragmentation 
function and the ALEPH data are shown in Fig. [?]. Furthermore, the results of this fit can 
be used to evaluate the photon + 1 jet rate for different values of y cut (y cut = 0.01, 0.1, 0.33). 
In each case, the calculated ^-distributions are compared with the leading order predictions 
and the ALEPH data in Fig. [?]. A good agreement is found for all values of y cut studied in the 
experimental measurement, reflecting the universality of the quark-to-photon fragmentation 
function. 

Since the present calculation is only lowest order in the strong coupling constant, there is 
a large ambiguity in the choice of a s and all values of a s are in principle equally valid. The 
value a s (Mz) = 0.124 used above was chosen such that the observed hadronic cross section 
can be reproduced by the 0(a s ) calculation. However, jet studies at LEP have indicated 
that lowest order calculations of jet observables can only be matched to experimental data 
if a larger value of the strong coupling constant to compensate for missing higher order 
contributions is used. In particular, the experimental data on the isolated photon + 2, 3 jet 
rates are well reproduced by a lowest order calculation M if a s (M z ) = 0.17. We therefore 
provide a fit of the quark-to-photon fragmentation function for this value of a s as well. We 
find 



^ ( r°Wo) = 




(10.3) 



with /i = 0.59 GeV. 

Figure [8] displays the fitted 0(aa s ) quark-to-photon fragmentation functions for a quark 
of unit charge in the MS-scheme at a factorization scale Hf — Mz> which is the only hard 
scale in the process. The corresponding 0(a) fragmentation function obtained in || is shown 
for comparison. It is apparent that these fragmentation functions differ largely for z — > 1, 
which indicates the need for a resummation of terms proportional to ln(l — z) in this region. 



10.2 The integrated rate for z > 0.95 

If the photon carries a fractional momentum z greater than 0.95 in the photon-jet it is 
considered to be isolated in the democratic clustering appoach used in the experimental 
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Figure 7: Comparison of the photon + 1 jet rate at leading and next-to-leading order with 
the ALEPH data. The non-perturbative quark-to-photon fragmentation function is fitted to 
the data for y cut = 0.06 only. The jet rates for the other values of y cut are then predictions 
from the leading order and next-to-leading order calculations. 



ALEPH analysis. In this approach the photon is treated like any other hadron when the 
clustering procedure is applied to form the jets. The division between isolated and non- 
isolated photons suggested by ALEPH in || is motivated by the fact that hadronization 
effects smear out the isolated photon peak from z — 1 to slightly smaller values of z. 

This definition of isolation is in contrast with that used in previous theoretical 0, 0] 
and experimental || analyses of isolated photon + n jet rates, where a two-step procedure 
was used to define an isolated photon. In these analyses, the photon was isolated from 
the other hadrons using a geometrical cone before these hadrons were clustered into jets. 
After the clustering had taken place the photon was required to be well-separated from all 
of the hadronic clusters and was said to be isolated if it was accompanied only by a small 
amount of hadronic energy inside the cone. In these theoretical calculations [|, ||, |7|], all of 
the quark-photon collinear and fragmentation contributions considered here were assumed 
to be negligible. As a result of these studies, it was found that large negative next-to-leading 
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Figure 8: The quark-to-photon fragmentation functions at leading and next-to-leading order 
as functions of z for a quark of unit charge. The factorization scale \xp is taken equal to Mz 
in both cases. 



order corrections were needed to provide a reasonable agreement between data and theory, 



precisely because soft gluons were excluded from the photon cone [11 



Using the fitted fragmentation function, we have evaluated the isolated photon + 1 jet 
rate in the democratic clustering approach up to 0(aa s ). The results of this calculation for 
the two different values of a s are compared with the data and the leading order prediction 
in Fig. H It appears that the leading-order curve provides an adequate description of the 
data and that the next-to-leading order curve improves the agreement between theory and 
experiment over the whole range of y cut - Furthermore, it is also apparent that the next-to- 
leading order corrections are moderate thereby demonstrating the perturbative stability of 
the definition of isolation within the democratic approach. 
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Figure 9: The integrated photon + 1 jet rate above z = 0.95 as function of y cut , compared 
with the full leading-order and next-to-leading order calculations including respectively the 
leading-order and next-to-leading order determined quark-to-photon fragmentation func- 
tions, D g _ >7 (z, fMp). 



11 Summary and Conclusions 

In this paper we have presented a complete 0(aa s ) calculation of the photon + 1 jet rate 
in e + e~ annihilation. Athough this calculation is formally only next-to-leading order in per- 
turbation theory, it contains several ingredients appropriate to next-to-next-to-leading order 
calculations of jet observables. For example, in addition to configurations where one final 
state particle is unresolved, there are also configurations where two particles are unresolved. 
In particular we encountered three double unresolved factors, corresponding to the triple 
collinear, soft/collinear and double single collinear limits of the 7* — > qqg'y subprocess. To 
analytically isolate the singularities associated with these configurations, we employed the 



resolved parton philosophy of |15| and introduced a theoretical parton resolution criterion 
y m i n . These previously unknown double unresolved factors are universal and are of the type 
expected to arise in any calculation of jet cross sections at next-to-next-to-leading order. 
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Matching the double unresolved parts of phase space with the single unresolved and resolved 
regions required some thought and was discussed at length in section ^. 

All these resolved and unresolved contributions from the double bremsstrahlung process 
must be combined with the single unresolved contributions from one-loop and fragmentation 
processes. Some single and double poles in e remain which have however exactly the right 
structure to be factorized in the next-to-leading order counterterm of the bare quark-to- 
photon fragmentation function. 

All remaining finite contributions can then be numerically evaluated for arbitrary jet 
clustering algorithms. The most stringent test on the consistency of our approach is provided 
by Fig. H where it is shown that the results of the numerical program are independent of the 
choice of the theoretical parameter y m ; n . This relies on the numerical cancellation of large 
logarithms of y min taken up to the third power. 

An important feature of this numerical program is the implementation of the hybrid 
subtraction method [S3] which is necessary to correctly treat the overlapping of two single 



collinear regions in the neighberhood of double unresolved regions of the four parton phase 
space. This overlapping was shown to be crucial for a correct and complete coverage of all 
singularities in the 7* — > qq^g process in section |3|. 

We have presented fits of the non-perturbative part of -D (J _ +7 for two values of a s (Mz) 
obtained by comparing our results with the ALEPH data f5|. In determining the fragmen- 
tation function, we have required that D q ^ is an exact solution of the 0(aa s ) evolution 
equation. This solution does not resum logarithms of the factorization scale, but does ensure 
that the photon + 1 jet rate is \ip independent. Furthermore, at large z, which corresponds 
to the most interesting 'isolated' photon part of the cross section, large logarithms of (1 — z) 
are present and should be resummed. 

As a first application of the fitted next-to-leading order quark-to-photon fragmentation 
function, we have compared the integrated 'isolated' photon rate for z > 0.95 with the 
same ALEPH data. The theoretical next-to-leading order predictions were found to describe 
the data well and to provide a better agreement between theory and experiment than the 
leading order calculation. Furthermore, the relatively small size of the next-to-leading order 
corrections indicates that the isolation definition used in this democratic clustering approach 
is perturbatively stable. 

Finally, the next-to-leading order quark-to-photon fragmentation function determined 
here is a universal function and appears in a variety of processes involving quarks and pho- 
tons. The most prominent examples of which are the prompt photon cross section at hadron 
colliders and the photon pair cross section at LHC. So far the fragmentation contributions 
to those processes have only been evaluated using model dependent assumptions for this 
fragmentation function |l7[ . It is reasonable to expect that the fragmentation function 
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determined directly from LEP data at high Q 2 and high z should significantly improve the 
theoretical predictions for these processes, and, in particular, may help to determine whether 
a Standard Model Higgs-boson of intermediate mass can really be detected via its two photon 
decay at the LHC. 
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